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cause of  its  excellent  neutron  economy,  its  versatility  with 
respect  to  fueling  strategies  and  the  high  fuel  burnup 
achievable  as  well  as  the  ability  to  approach  continuous 
fueling  operation,  the  'pebble  bed  HTGR  was  judged  to  be  the 
best  suited,  of  reactors  currently  of  commercial  interest, 
for  implementation  of  the  breed/burn  fuel  cycle  in  a  thermal/ 
epi thermal  system. 

Three  necessary  (but  not  sufficient)  criteria  for  success- 
ful system  operation  were  defined:   (1)  achieving  system  criti- 
cality,  k^  >^  1.0;  (2)  achieving  a  breed  zone  fissile  production 
rate  such  that  the  asymptotic  enrichment  attained  in  the  breed 
zone  and  the  end-of-life  burnup  in  the  burn  zone  (expressed 
in  percent)  are  related  by  Cp^^Y   1  ^EOL'  ^^^  ^-^^    accounting  for 
zone-to-zone  leakage  by  providing  adequate  excess  reactivity, 
shown  to  require  that  k   >  1.5. 
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State-of-the-art  computation  methods,  PDQ -7/A2'^PX,  were 
employed  for  the  analysis  of  homogeneous  breed/burn  pebble  bed 
}ITGR  regions.   Breed  zones  achieved  their  highest  asymptotic 
enrichment  ('v-3.5%)  for  carbon-to-heavy  metal  atom  ratios, 
C/IIM  =  110;  when  charged  to  the  burn  zone  this  maximized 
koo  (^^1.3).   Highly-m.oderatcd  burn  zones  (C/HM  =;  450)  were  most 
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were  able  to  achieve  the  k^^  (k„  >^  1.5)  required  to  overcome  the 
effect.j  of  zone- to- zone  leakage.   For  these  compositions  then, 

successful  operation  of  a  breed/burn  system  was  not  possible. 
Further  study  of  reactor  systems  operating  with  very  hard  breed 
zone  eipthermal  spectra  (C/HM  <  110)  should  be  undertaken  if  a 
feasible  breed/burn  system  is  to  be  pursued  since  it  was  found 
that  asymptotic  zone  enrichment  reaches  a  maximum  in  a  spectrum 
between  that  typical  of  a  hard-spectrum  HTGR  and  soft-spectrum 
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ABSTRACT 


The  feasibility  of  the  breed/burn  fuel  cycle  concept, 
originally  proposed  for  Pu239/U  FBR  applications,  has  been 
evaluated  for  use  in  a  U233/Th  fueled  pebble  bed  HTGR.   Be- 
cause of  its  excellent  neutron  economy,  its  versatility  with 
respect  to  fueling  strategies  and  the  high  fuel  burnup 
achievable  as  well  as  the  ability  to  approach  continuous 
fueling  operation,  the  pebble  bed  HTGR  v/as  judged  to  be  the 
best  suited,  of  reactors  currently  of  commercial  interest, 
for  implementation  of  the  breed/burn  fuel  cycle  in  a  thermal/ 
epithermal  system. 

Three  necessary  (but  not  sufficient)  criteria  for  success- 
ful system  operation  were  defined:   (1)  achieving  system  criti- 
cality,  k   >^  1 . 0 ;  (2)  achieving  a  breed  zone  fissile  production 
rate  such  that  the  asymptotic  enrichment  attained  in  the  breed 
zone  and  the  end-of-life  burnup  in  the  burn  zone  (expressed 
in  percent)  are  related  by  e^^sY  ^  ^EOL'  ^^"^^  (^)  accounting  for 
zone-to-zone  leakage  by  providing  adequate  excess  reactivity, 
shown  to  require  that  k   >  1.5. 

state-of-the-art  computation  methods,  PDQ  -7/AI'lPX,  were 
employed  for  the  analysis  of  homogeneous  breed/burn  pebble  bed 
HTGR  regions.   Breed  zones  achieved  their  highest  asymptotic 
enrichment  ('^>3.5%)  for  carbon-to-heavy  metal  atom  ratios, 
C/HM  =  110;  when  charged  to  the  burn  zone  this  maximized 
koo  (^^1.3).   Highly-moderated  burn  zones  (C/HM  ^    450)  were  most 
able   to  attain  high  fuel  burnups  (^^100,000  MWD/MTHM)  while 
sustaining  the  highest  possible  values  of  koo  .   No  combination 
of  the  compositions  evaluated  here  (110  <_   C/HM  >_   450)  ,  however, 
were  able  to  achieve  the  k^  (k^o  >^  1.5)  required  to  overcom.e  the 
effects  of  zone-to-zone  leakage.   For  these  compositions  then. 
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successful  operation  of  a  breed/burn  system  was  not  possible. 
Further  study  of  reactor  systems  operating  with  very  hard  breed 
zone  eipthermal  spectra  (C/HM  <  110)  should  be  undertaken  if  a 
feasible  breed/burn  system  is  to  be  pursued  since  it  was  found 
that  asymptotic  zone  enrichment  reaches  a  maximum  in  a  spectrum 
between  that  typical  of  a  hard-spectrum  HTGR  and  soft-spectrum 
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I.   INTRODUCTION 

1 . 1   Foreword 

In  the  past  the  general  public  has  understood  the  term 
"nuclear  proliferation"  in  a  restricted  sense,  referring  almost 
exclusively  to  the  widespread  deployment  of  nuclear  weapons 
throughout  the  world.   The  realization  that  a  growing  number  of 
foreign  countries  are  operating  commercial  nuclear  powered 
electrical  generating  stations  and  the  implied  relationship  be- 
tween nuclear  electric  plants  and  weapons  has  resulted  in  a 
broader  interpretation  of  proliferation.   In  the  United  States, 
this  interpretation  has  achieved  the  status  of  official  govern- 
ment policy  and  has  led  to  both  national  and  international 
reassessments,  e.g.  the  Nonprolif eration  Alternative  System 
Assessment  Program  (NASAP)  and  the  International  Fuel  Cycle 
Evaluation  (INFCE)  efforts.   The  question  of  the  safeguarding 
of  sensitive  nuclear  materials  has  become  a  major  area  of 
emphasis  in  furthering  the  development  and  application  of  the 
nuclear  fuel  cycle  for  commercial  electrical  production. 

Previously  it  was  envisioned  that  light  water  reactors 
(LVJR)  would  serve  an  interim  role  as  precursors  to  fast  breeder 
reactors  (FBR) .   This  status  is  testified  to  by  the  omission  in 
most  of  the  older  long  range  systems  studies  of  "advanced  LWRs . " 
Until  such  time  as  FBRs  could  become  commercially  feasible  and 
sufficient  fissile  reserves  could  be  bred  to  permit  their 
widespread  deployment,  LWRs  were  to  be  the  mainstay  of  commer- 
cial nuclear  power. 
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Currently,  however,  the  transition  to  a  breeder  economy 
is  not  envisioned  for  the  United  States,  at  least  in  the 
short  term,  so  it  is  prudent  to  consider  other  options  to 
meet  the  needs  for  a  growing  electrical  generating  capacity 
[Nl] .   Since  the  objectionable  feature  of  the  FBR  is  that  it 
utilizes  fertile  uranium-238  (U238)  to  breed  large  amounts  of 
fissile  Plutonium,  new  options  must  permit  utilizing  nuclear 
fuels  in  such  a  manner  that  the  production  of  and/or  commer- 
cial traffic  in  materials  suitable  for  weapons  applications  is 
minimized  at  all  levels. 

The  dilemma  posed  by  this  objective  is  highlighted  by 
rising  oil  prices  and  limited  and  unevenly  distributed  national 
energy  resources.   Many  nonnuclear  weapons  states,  desiring  to 
provide  for  their  own  growing  energy  needs,  have  opted  to 
develop  civilian  nuclear  power  programs.   In  the  United  States, 
the  lack  of  expansion  of  enrichment  capacity  or  the  commercial 
deployment  of  reprocessing  and  mixed  oxide  fabrication  capabil- 
ity has  increased  the  economic  and  commercial  incentives  for 
other  countries  to  develop  such  capabilities  independently  [PI] . 

In  view  of  the  above  considerations,  development  and  evalu- 
ation of  the  commercial  feasibility  of  schemes  utilizing  the 
nuclear  fuel  cycle  in  a  manner  so  as  to  minimize  proliferation 
risks,  while  providing  suitable  economic  and  commercial  incen- 
tives for  their  widespread  deployment,  has  become  an  important 
national  and  international  objective.   It  is  the  purpose  of  the 
present  thesis  to  evaluate  one  particular  class  of  prolifera- 
tion-resistent  fueling  schemes;  here  designated  the  "breed/burn" 


-13- 


concept.   As  will  be  discussed  further,  in  its  fullest  develop- 
ment this  scheme  would  allow  charging  only  fertile  material  to 
a  reactor  after  initial  startup  and  dispatch  of  fully  spent 
fuel  to  disposal  rather  than  reprocessing  and  reuse. 
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1 .  2   Background 

Throughout  its  history,  commercial  nuclear  power  in  the 
United  States  has  been  subjected  to  several  periods  of  chang- 
ing emphasis.   To  the  public  utilities,  short  term  fuel  cycle 
costs  have  always  been  a  most  important  consideration.   While 
long-term  trends  are  reviewed  and  evaluated,  the  utility 
always  seeks  to  minimize  the  current  costs  of  electricity  to 
the  consumer,  hence  causing  the  short-term  economic  consider- 
ations to  dominate.   Indeed,  the  very  practice  of  present 
value  discounting  emphasizes  the  short-term. 

National  energy  planning,  on  the  other  hand,  is  concerned 
with  being  able  to  insure  the  availability  of  adequate  energy 
resources  within  reasonable  economic  limits  for  the  very  long- 
term.   With  specific  application  to  the  nuclear  fuel  cycle, 
although  most  quantities  of  high-cost  fissile  uranium  in  the 
earth's  crust  can  conceivably  be  considered  a  source  of  energy, 
the  cost  of  utilizing  it  is  prohibitive.   In  this  light,  a 
more  economically  viable  solution  to  the  long-term  energy 
supply  problem  is  to  produce  most  future  fissile  materials  by 
irradiating  fertile  materials  in  the  fuel  elements  of  a  nuclear 
reactor.   However,  to  accomplish  this  without  relying  on  the 
commercialization  of  the  FBR,  it  becomes  necessary  to  operate 
other  types  of  reactors  as  near-  or  net-breeders. 

Current  pressurized  water  reactors  (PWR)  operate  with 
conversion  ratios  on  the  order  of  0.60  [Dl]  on  the  uranium-238 
fuel  cycle.   High  temperature  gas-cooled  reactors  (HTGR) ,  how- 
ever, typically  operate  with  thorium-232  (Th232)  to  uranium-233 
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(U233)  conversion  ratios  of  0.76  [Dl],  based  on  optimization 
of  near  term  economics.   Conversion  ratios  of  0.84  are  feasible 
and  favored  for  ore  conservation,  however  [Dl] ,  and  even 
higher  values  are  readily  obtainable  [L4,R2,S3].   Thus,  we  will 
focus  on  this  inherently  superior  class  of  reactors  in  the  pre- 
sent study.   Most  HTGRs  built  or  designed  to  date  utilize  a 
highly  enriched  uranium/ thorium  (HEU/Th)  fuel  cycle  which 
necessitates  uranium-235  (U235)  enrichments  as  high  as  93%. 
Obviously,  to  meet  nonprolif eration  objectives,  alternatives 
to  this  fuel  cycle  must  be  evaluated. 

The  HTGR  coated-particle  fuel  design  allows  a  variety  of 
alternative  cycles,  each  having  characteristic  advantages  and 
disadvantages  in  particular  circumstances.   Designing  alterna- 
tive fuel  cycles  that  can  use  the  same  core  layout,  fuel  zoning 
and  fuel  blocks  as  the  HEU/Th  cycle  means  that  conversion  from 
one  fuel  cycle  to  another  can  be  accomplished  in  available 
plant  designs  merely  by  replacing  fuel  blocks  in  each  segment 
with  alternative  cycle  blocks  during  the  course  of  the  regular 
refueling  sequence.   Since  the  parameters  of  fissile-to-fertile 
and  heavy  metal-to-moderator  atom  ratios  control  the  conversion 
ratio  of  the  reactor  system,  the  versatility  of  the  HTGR  makes 
it  a  logical  choice  for  tailoring  fuel  loadings  and  burnups , 
based  on  changing  economic  and  resource  availability  conditions, 

Even  more  versatile  and  convenient  for  implementation  of 
alternative  fuel  cycles  among  the  advanced  HTGRs  is  the  pebble 
bed  reactor  [L4,Tl].   The  injection-molded  spherical  graphite 
fuel  balls  used  to  fuel  this  system  can  easily  be  fabricated 
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with  any  desired  fissile  or  fertile  fuel  density,  and  immedi- 
ately be  introduced  into  the  continuous  refueling  cycle. 

Operating  an  HTGR  as  a  near-to-net-breeder  utilizing  the 
thorium  fuel  cycle  requires  low  power  density  operation  with 
heavy  thorium  loadings.   Since  operating  trends  are  dictated 
by  short  term  energy  costs,  obtaining  a  maximum  conversion 
ratio  and  better  fuel  utilization  is  not  necessarily  a  desirable 
near  term  goal.   In  obtaining  a  higher  conversion  ratio,  in- 
creasing the  thorium  loading  in  the  core  while  decreasing  the 
core  power  density  and  fuel  residence  time  would  incur  off- 
setting economic  penalties  which  in  turn  may  or  may  not  be 
balanced  by  the  advantages  of  improved  fuel  utilization.   The 
penalties  would  depend  on  projected  economic  conditions,  the 
projected  availability  of  uranium  ore,  costs  for  separative 
work,  reprocessing  and  recycle,  and  transportation  and  disposal 
costs. 

With  increasing  thorium  loading,  although  an  increase  in 
the  conversion  ratio  will  be  realized,  a  higher  fissile  inven- 
tory is  also  required.   The  resulting  higher  initial  costs  will 
be  attractive  if  sufficient  money  is  saved  over  the  life  of 
the  reactor.   This  outcome  is  favored  in  an  economic  environment 
where  large  price  increases  are  projected  for  fresh  feed  fuel. 

Lowering  the  core  power  density  in  an  HTGR  will  result  in 
a  larger  core  volume  for  a  given  thermal  output.   Hence,  addi- 
tional volume  is  available  for  fertile  thorium  loading,  pro- 
ducing the  same  net  result  as  previously  described.   Addition- 
ally, neutron  utilization  is  improved  due  to  a  reduction  in 
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losses  to  xenon-135  (Xel35)  and  protactinium-233  (Pa233)  [K2]. 
The  larger  reactor  volume  will  increase  initial  capital  invest- 
ment and  will  also  change  the  required  coolant  pumping  power 
of  the  system.   Keeping  the  coolant  temperature  increase  con- 
stant across  the  core  while  maintaining  similar  height-to- 
diameter  ratios,  the  larger,  low  power  density  core  will  require 
a  higher  coolant  mass  flow  for  a  given  thermal  output.   Because 
the  HTGR  uses  helium  as  its  primary  cooling  medium,  very  high 
coolant  mass  flow  rates  require  that  a  significant  fraction  of 
the  plant  power  output  be  utilized  in  operating  the  coolant 
pumps.   Therefore,  a  variation  in  the  core  power  density  that 
will  increase  the  pumping  requirements  can  add  a  sizeable 
financial  burden  to  the  overall  operating  costs. 

In  addition  to  the  HEU/Th  fuel  cycle,  the  use  of  a  uranium- 
233/thorium  (U233/Th)  cycle,  a  lov;-enriched  uranium  (LEU)  cycle 
or  various  plutonium/uranium  (Pu/U)  or  plutonium/uranium/ thorium 
(Pu/U/Th)  fuel  cycles  are  also  possible  in  the  KTGR  [Bl,S3,Tl]. 
For  all  of  the  options  the  reduced  total  uranium  ore  requirement 
achieved  by  increasing  the  conversion  ratio  requires  accepting 
either  a  higher  capital  investment  or  higher  operating  costs 
in  other  areas . 

The  present  work  focuses  primarily  on  the  U233/Th  cycle 
(because  of  its  well  known  superiority  in  thermal  and  epi- 
thermal  systems)  and  the  possibility  of  its  utilization  in  a 
system  designed  to  accomplish  most  key  nonprolif eration 
objectives.   Since  U233  is  required  as  the  initial  fissile 
loading  in  this  cycle,  the  resource  conservation  is  not  direct. 
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but  rather  indirect.   To  provide  the  initial  U233  loading, 
adequate  fissile  production  in  prebreeders,  operating  expressly 
for  the  production  of  U233  is  required.   To  accomplish  the 
objectives  of  nonprolif eration  effectiveness  we  will  evaluate 
the  feasibility  of  using  the  U233/Th  cycle  in  either  the 
internally  recycled  or  the  "Once-Through-Then-Out"  (OTTO) 
modes  in  a  pebble  bed  reactor  system  [Tl] . 

In  the  projection  of  any  long-term  energy  strategy,  the 
extension  of  uranium  resources  must  be  considered  as  most 
important,  even  in  view  of  the  increasing  relevance  of  the 
objectives  of  nonprolif eration.   The  present  work,  however, 
will  attempt  only  to  evaluate  the  feasibility  of  a  single 
scheme  for  effectively  providing  a  viable  solution  to  the 
proliferation  problem  of  safeguarding  sensitive  nuclear 
reactor  materials  through  all  steps  of  the  fuel  cycle  without 
regard  to  its  long-term  effects  on  either  the  utilization  of 
resources  or  the  economics  of  implementing  such  a  scheme 
commercially.   However,  because  breakeven  breeding  is  neces- 
sary to  fully  achieve  the  nonprolif eration  objectives,  long- 
term  ore  utilization  should  be  quite  efficient. 

1.3   Previous  Work 

The  use  of  the  breed/burn  concept  was  proposed  by  Fischer 
et  al  [Fl]  for  implementation  in  a  reactor  operating  with  a 
hard  spectrum  core.   Depleted  uranium  monocarbide  feed  fuel  would 
be  introduced  to  breed  fissile  plutonium,  and  after  a  period 
of  residence  in  the  hard  spectrum  the  fuel,  now  slightly 
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enriched,  v/ould  be  moved  to  the  more  thermal! zed  outer  region  of 
the  fast  core.   Energy  would  be  generated  as  the  fuel  burns 
to  a  design-dependent  burnup  limit  while  the  fissile- to-fertile 
fuel  ratio  is  reduced  toward  the  equilibrium  value  for  that 
neutron  energy  spectrum. 

The  present  work  extends  the  concept  of  the  breed/burn 
cycle  to  assess  the  feasibility  of  its  implementation  in  the 
thermal/epi thermal  spectrum  of  the  HTGR.   Nonprolif eration 
objectives  of  reducing  fissile  enrichments  can  be  realized 
in  addition  to  the  benefits  of  decoupling  this  type  of  reactor 
from  breeders  and  the  constraints  related  to  reprocessing 
facilities.   The  exception,  of  course,  is  the  need  to  obtain 
the  initial  startup  core,  here  requiring  a  substantial  U233 
inventory. 

The  introduction  of  this  fueling  scheme  into  the  pebble 
bed  reactor  system  provides  the  ability  to  readily  adjust  the 
core  neutron  spectrum  simply  by  varying  the  moderator-to- 
fuel  (carbon-to-heavy  metal)  ratio.   Additionally,  core 
zoning  (establishing  separate  breed  and  burn  regions)  can 
easily  be  accomplished  in  this  application.   A  given  ball, 
during  its  lifetime,  would  reside  in  two  distinct  regions. 
Radial  zoning  would  make  coolant  temperature  control  difficult. 
Since  the  breed  zone  would  operate  at  a  lower  power  density 
than  the  burn  zone,  cold  coolant  bypass  v/ould  degrade  the  bulk 
core  exit  temperature.   Therefore,  the  core  will  be  divided 
into  axial  zones,  allowing  it  to  be  seen  that  for  a  fertile  load 
ball,  the  "breed  zone"  can  be  defined  as  that  region  of  the 


reactor  in  which  the  fissile  concentration  of  the  ball  is  in- 
creasing.  Similarly,  the  "burn  zone"  is  that  region  in  which 
the  fissile  concentration  decreases.   The  desired  breeding  and 
burning  effects  can  be  accomplished  by  operating  the  breed  zone 
in  a  predominantly  epithermal  spectrum.   Thus,  in  its  lifetime, 
a  fertile-fueled  ball  would  travel  through  both  zones,  enter- 
ing the  thermal  burn  zone  (together  with  pure  graphite  balls 
necessary  to  attain  the  desired  moderator-to-fuel  ratio) 
after  first  passing  through  the  breed  zone  (which  contains 
only  fueled  balls) .   A  schematic  diagram  of  system  operation 
is  illustrated  in  Figure  [1.1]. 

In  successive  passes  through  the  breed  zone,  the  concentra- 
tion of  U233  is  built  up  to  the  maximum  possible  amount  for 
the  given  spectrum.   It  is  then  cycled  (with  additional  gra- 
phite-only balls)  into  the  burn  zone,  again  being  recycled 
for  successive  passes  until  a  design-limited  burnup  is 
attained.   At  this  point  the  depleted  balls  are  removed  from 
the  system  for  disposal. 

Operating  on  this  type  of  fuel  cycle,  the  proposed  reactor 
would  allow  very  high  burnups,  and  as  a  result,  the  total 
fissile  enrichment  of  the  discharged  balls  would  be  very  low, 
particularly  if  the  number  of  passes  could  be  minimized, 
in  which  case  discharge  compositions  approaching  those  of  the 
OTTO  cycle  could  be  envisioned.   If  the  reactor  core  can  be 
adequately  segregated  so  as  to  maintain  distinct  breed  and 
burn  regions,  then  the  practicality  of  achieving  the  proper- 
ties sought  by  Fischer  et  al  can  be  evaluated. 
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Work  on  HTGRs  has  been  devoted  mainly  to  the  use  of  the 
HEU/Th  fuel  cycle,  both  in  traditional  graphite  block-element 
reactors  and  pebble  bed  systems.   Many  alternative  cycles, 
previously  noted,  have  also  been  investigated  more  recently 
and  in  less  detail,  both  from  their  ability  to  meet  nonpro- 
liferation  objectives  as  well  as  their  respective  roles  in  the 
extension  of  uranium  resources.   Current  studies  [D2,T1,V1] 
are  focusing  on  the  flexibility  of  the  HTGR  fuel  cycle,  which 
permits  it  to  adapt  to  variations  in  fuel  materials  relatively 
efficiently.   While  favorable  in  many  respects,  this  feature 
also  makes  it  difficult  to  determine  optimal  strategies.   This 
extra  flexibility,  relative  to  other  reactor  types,  results 
from  the  particle  loading  system,  which  allows  one  to  vary 
the  fertile  and  fissile  loadings  almost  independently  and  also 
to  adjust  the  resonance  shielding,  not  only  by  varying  the 
size  and  loading  of  the  fuel  elements,  but  also  by  adjusting 
the  dimensions  of  the  individual  fuel  kernels  [Dl] . 

Utilizing  this  flexibility  it  is  possible  to  operate 
HTGRs  with  conversion  ratios  slightly  greater  than  unity 
[R2,S3].  But  considering  resources  and  economics,  the  high 
fissile  inventories  required  and  the  resulting  recycle  losses 
have  made  fuel  designs  with  conversion  ratios  on  the  order 
of  0.9  appear  more  promising  [S3],  because  near  this  value, 
burnup  and  processing  losses  reach  a  minimum  [Dl] . 

Fuel  cycle  optimization  in  the  HTGR  has  traditionally  been 
most  heavily  concerned  with  minimizing  fuel  cycle  cost.   It  is 
clear  that  this  approach  does  not  fully  utilize  the  potential 
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in  the  HTGR  for  resource  conservation  through  the  use  of  the 
U233/Th  fuel  cycle.   To  obtain  a  breeding  system  with  this 
cycle,  if  less  than  r]-2   neutrons  are  lost  to  parasitic  ab- 
sorbers, then  after  one  neutron  is  utilized  in  sustaining 
criticality,  one  additional  neutron  is  available  for  absorp- 
tion by  the  fertile  thorium  to  produce,  after  3-decay,  a 
new  fissile  U233  atom.   The  selection  of  the  U233/Th  cycle 
for  this   application  is  clear  when  the  spectrum-averaged 
neutron  yield  per  absorption,  shown  in  Table  [1.1],  is  con- 
sidered for  this  reactor  type.   For  various  carbon-to-heavy  metal 
atom  ratios  (N  /N^^j^)  only  fissile  U233  has  a  value  of  n  ade- 
quately high  to  allow  breeding. 

The  high  flexibility  in  the  design  of  fuel  elements  and 
the  reactor  refueling  scheme  in  the  pebble  bed  HTGR  has  made 
its  use  as  a  thermal  breeder  possible  [R2] .   Using  highly 
enriched  U233  ( (U233/U) -  0  . 9 )   produced  in  a  prebreeder  as  feed 
fuel,  thick  radial  blankets  of  thorium,  high  heavy  metal  load- 
ings in  the  fuel  elements  and  relatively  low  heavy  metal  burnups 
will  allow  the  pebble  bed  system  to  be  self-supporting;  inde- 
pendent of  uranium  resources  after  the  first  core  inventories 
have  been  built  up.   Such  designs  have  been  shown  to  achieve 
conversion  ratios  greater  than  1.05  [R2] . 

Other  studies  [T1,V1],  dealing  more  with  nonprolif eration 
and  ore  conservation  strategies  than  with  breeding  potential 
have  shown  that  proliferation  resistance  is  best  achieved  in 
the  lov7-enriched  mixed  oxide  cycle  (mixed  uranium/ thorium 
oxides  with  U235/U-  0  .  20)  .   Reactor  loading  is  restricted  to 
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TABLE  1.1 
SPECTRUM-AVERAGED  NEUTRON  YIELD 
PER  FISSILE  ABSORPTION* 

CARBON-TO-HEAVY  METAL 
FISSILE        ATOM  RATIO  (N^^/N   )  + 
ISOTOPE      110     250     325    450 
U233       2.209   2.240   2.246   2.253 
U235       1.928   1.989   2.001   2.011 
Pu239  ■     1.775   1.788   1.791   1.794 


va 


f 


a 
a 


"^  C/HM  ratios  correspond  to  spectrum  averaged  cross- 
section  data  given  in  Appendices  B  and  D. 
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non-sensitive  materials  while  95%  of  the  bred  fissile  plutonium 
is  burned  in  situ  [Tl] . 

Resource  conservation,  however,  is  most  favorable  for  OTTO 
cycles  using  highly-enriched  U/Th  mixed  oxide  (  (U235/U)  -0.93) 
fuels  designed  for  burnups  of  100  to  130  GWD/MTHM  [Tl] . 

From  the  above  discussion,  it  can  be  seen  that  the  wide 
latitude  of  compositions  and  fueling  schemes  available  in  the 
HTGR  allow  considerable  breadth  in  designing  strategies  empha- 
sizing objectives  of  proliferation  resistance,  fissile  breeding 
potential,  uranium  resource  conservation  or  fuel  cycle  econom- 
ics.  Each  strategy  has  characteristic  advantages  and  disadvan- 
tages, depending  on  the  design  criteria  for  the  chosen  fuel 
cycle.   The  breed/burn  concept,  when  applied  to  the  pebble  bed 
HTGR  system,  combines  the  objectives  of  nonprolif eration  and 
fissile  breeding.   It  is  anticipated  that  the  HTGR  system  can 
best  provide  the  versatility  needed  to  devise  such  a  scheme  for 
a  thermal/epithermal  system. 

1 . 4   Outline  of  Present  Work 

The  main  body  of  this  report  consists  of  three  chapters. 
Chapter  2  discusses  the  operation  of  the  breed/burn  fuel  cycle 
in  the  pebble  bed  system  and  the  relationship  of  the  neutronic 
parameters  of  the  two  regions.   Additionally,  an  evaluation  is 
made  of  the  potential  of  the  pebble  bed  system  for  successful 
operation  utilizing  the  breed/burn  concept  using  one  group 
models.   Discussion  is  presented  of  the  effects  of  spectrum 
hardening  on  the  rate  and  amount  of  fissile  breeding  obtainable 
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Chapter  3  describes  the  computational  methods  used  for 
evaluation  of  the  reactor  system,  nuclear  data  input,  cross- 
section  generation  and  the  reactor  model. 

The  results  of  the  analysis  are  discussed  in  Chapter  4 
with  regard  to  reactor  criticality,  burnup  limits  and  rates 
of  fissile  production  as  a  function  of  composition.   Possible 
modifications  to  the  system  are  suggested. 

Chapter  5  sumniarizes  the  analysis  and  reiterates  the  main 
conclusions.   Recoirimendations  are  made  for  future  work.   Appen- 
dices containing  supplementary  information  are  also  included. 
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II.   BREED/BURN  OPERATION 

2 . 1  Introduction 
The  pebble  bed  HTGR  is  expected  to  provide  the  best  avail- 
able system  for  implementation  of  the  breed/burn  concept  into 

a  thermal/epithermal  reactor  because  of  it  compositional 
flexibility  and  high  burnup.   In  order  to  utilize  the  breed/ 
burn  fuel  cycle  in  this  system,  the  feasibility  of  operating  it 
as  a  two-zone  reactor  will  be  investigate.   Dividing  the 

core  axially  will  allow  the  two  zones  to  exist  in  harmonv  with 

I 

'  the  flow-through  fuel  cycle,  such  that  all  fuel  elements  will 

move  through  both  zones  during  successive  passes  through  the 

core. 

I 

The  purpose  of  this  chapter  is  to  consider  the  effects  of 

various  parameters  on  the  rate  of  U233  production  in  the  breed 
zone  of  such  a  system  and  to  evaluate  the  potential  for  success- 
ful operation  in  the  pebble  bed  HTGR.   Possible  modifications 
will  be  discussed  to  enhance  system  performance. 

2 . 2  Breed/Burn  Potential  for  U233  Production 

The  objective  is  to  determine  if,  in  fact,  it  is  a  feasible 
consideration  to  attempt  using  the  breed/burn  system  in  a  peb- 
ble bed  HTGR.   For  this  analysis,  the  system  considered  will 
alternate  breed  and  burn  regions  axially  for  the  length 
of  the  reactor.   The  pebble  bed  HTGR  is  a  very  large,  low 
power  density  system,  and  hence  the  net  neutron  leakage  from 
the  core  is  relatively  low.   By  varying  the  composition 
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and  geometry  of  each  region,   the  respective  neutrohic  charac- 
teristics of  the  zones  can  !oe   altered,  thereby  changing  the 
internal  leakage  between  zones  and  hence  changing  the  operating 
envelope  of  the  entire  system. 

First,  an  analytical  method  will  be  used  to  evaluate  the 
parameters  having  the  greatest  effect  on  the  rate  of  fissile 
breeding,  and  then  the  entire  system  will  be  evaluated  for  its 
ooerational  capabilities. 

2 . 3   Factors  Affecting  the  Breeding  Ratio 

The  customary  index  of  fissile  breeding  (here  of  U233 
from  neutron  capture  in  fertile  Th232)  is  represented  by  the 
breeding  ratio,  b,  defined  as: 


_  Rate  of  fissile  production  in  the  system. 
Rate  of  fissile  consumption  in  the  system 


or. 


■     C  °2  +  C   02 

A     +  A, 
c      bz 


where       C  =  neutron  capture  rate  in  the  species  indicated, 

A  =  neutron  absorption  rate  in  the  species 
indicated, 

c  =  burn  zone 


and         bz  =  breed  zone. 

For  simplicity  in  this  analysis,  the  effects  of  only  the 
most  significant  materials  will  be  considered.   While  it  is 
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realized  that  a  rigorous  approach  which  would  include  the  pro- 
duction and  depletion  of  protactinium  and  U235  and  the  deple- 
tion of  Th232   could  yield  more  precise  conclusions,  the 
factors  controlling  the  rate  of  fissile  U233  production  can 
be  examined  adequately  by  this  more  simple  analysis. 

In  each  region  a  neutron  balance  can  be  written  as: 

Production  =  Losses  +  Absorption  +  Leakage 
which,  for  the  burn  zone,  is  given  by: 


VF  "  +  ^p  02  ^  ^  ^p  23^p  02^^  23^^  02  +  A  ^         (2.2) 
c      C       cccccc 


similarly  for  the  breed  zone: 


VF   "  +  VF   02  ^  L   =  F,  "+F   °2  +  c   23^C   02+.   P  ,2.3, 

bz       bz       c     bz      bz      bz      bz      bz 

where 

F  =  fission  rate  in  the  species  indicated, 

L  =  neutron  leakage  from  the  region  indicated 

P 
and  A   =  parasitic  absorption. 

Combining  these  terms  the  breeding  ratio  of  the  system  can  be 
expressed  as: 

(V-I)  F  02    A^ 

^S  S 
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where   n^  -  mean  neutron  yield  per  fissile  absorption  in  the 
system 


23 
F 
s 

A  23 
s 


(2.5) 


and      s  =  sum  of  indicated  interactions  in  the  breed 
and  burn  zones. 

Similarly,   the  breeding  ratios ' for  the  burn  zone  alone 
and  the  breed  zone  alone,  respectively,  are  given  by: 

(v-l)  F  °2    ;,  P 

(v-l)  F^  02       P 

^"<5      \z    ~-    %z  *   ,  23      -  rfe  -  1  (2.7) 

A  A 

s  s 

—    23 

where       n^  =   ^j-  (2.8) 

A 
s 

-  r.     23 

"^  ^bz 
and      n^^  =    ^23   "    .  ^2. 9) 


Equations  (2.4),  (2.6)  and  (2.7)  are  composed  of  three  main 
terms.   The  first  term,  t]  ,    dominates  the  expressions,  and 
hence  increasing  n  will  increase  the  breeding  ratio  proportion- 
ately.  This  can  be  accomplished  by  hardening  the  core  neutron 
spectrum  through  an  increase  in  the  concentration  of  the 
fertile  and  fissile  species.   That  is,  lowering  the  carbon- to- 
heavy  metal  (C/HM)  atom  ratio  and  thus  decreasing  the  modera- 
tion results  in  a  harder  neutron  spectrum. 
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The  second  term,  (v-1)  F   /A^"^,  represents  the  gain  in 
n  due  to  fast  fissioning  of  fertile  material.   Because  of  the 
small  threshold  fission  cross-section  of  Th232  this  term  is 

relatively  small. 

P   23 
The  third  term,  A  /A   ,  represents  the  neutrons  lost  to 

parasitic  absorption  and  leakage  and  hence  is  negative.   Thus, 
the  objective  of  increasing  the  amount  of  fissile  breeding 
necessitates  minimizing  the  amount  of  parasitic  absorption 
and  leakage  losses.   Fuel  material  absorption  cross-sections 
and  those  of  structural  materials,  therefore,  become  important 
design  considerations.   From  this  analysis  we  must  conclude 
that  in  order  to  increase  the  breeding  ratio  in  our  system   a 
hardened  neutron  spectrum  is  generally  favored  since  it  mini- 
mizes parasitic  absorption  and  maximizes  neutron  yield,  as  will 
be  shown  in  succeeding  sections.   Leakage  losses  may  or  may 
not  be  improved,  however.   The  first  design  strategy  clue  is 
then  to  obtain  as  hard  a  spectrum  as  practical  in  both  the 
breed  and  burn  zones  within  other  constraints.   Because  we 
want  to  be  able  to  produce  adequate  masses  of  fissile  material 
in  the  breed  zone  to  fuel  the  burn  zone,  our  attention 
focuses  on  the  parameters  of  the  breed  zone.   It  is  realized 
at  this  point,  however,  that  the  two  zones  will  have  significant 
coupling  in  their  characteristics.   That  is,  considering  a 
given  number  of  neutrons  in  the  system,  favoring  neutron  avail- 
ability in  one  region  will  decrease  the  effectiveness  of  the 
other.   If  the  breed  zone  is,  under  any  conditions,  unable  to 
breed  adequate  masses  of  U233,  then  the  neutron  balance  in  the 
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burn  zone  will  still  be  unable  to  effect  satisfactory  system 
operation.   It  is  for  this  reason  that  the  breed  zone  para- 
meters will  be  examined  in  more  detail. 

2.4   Breed  Zone  U233  Conversion 

To  some  degree,  breeding  of  U233  will  occur  throughout 
the  reactor.   In  applying  the  breed/burn  concept  to  the  pebble 
bed  HTGR  a  criterion  for  successful  operation  has  been  set 
that  adequate  breeding  in  the  breed  zone  alone  can  be  achieved 
to  provide  sufficient  total  fissile  mass  to  fuel  the  burn  zone. 
Obviously,  a  simultaneous  objective  is  to  maintain  system 
criticality  for   a  reasonable  cycle  lifetime.   Breeding  of  U233 
in  the  entire  system  can  be  divided  into  two  parts,  correspond- 
ing to  burn  zone  breeding  and  breeding  in  the  breed  zone. 

For  the  burn  zone: 


,   _  Rate  of  fissile  production  in  the  burn  zone - 

c    Rate  of  fissile  consumption  in  the  breed  and  burn  zones 

C  02 
^  (2.10) 


23       23  ' 
(A  "^^    +  A,  ^^) 
c  bz   ' 


while  for  the  breed  zone; 


_  Rate  of  fissile  production  in  the  breed  zone 

bz  ~  Rate  of  fissile  consumption  in  the  breed  and  burn  zones 

'^^  .  (2.11) 
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Solving  Eq .  (2.3)  for  the  neutron  capture  rate  in  Th233 
in  the  breed  zone: 


r       02  =  vF^  23  +  vF^  02  -  F^  23      02      23      P 
bz        bz       bz      bz      bz      bz      bz     c 


=  (V  -  1)  (F,   ^^  +  F,  °^)  -  a  ^"^  -  A^  ^  +  L 
bz      bz       bz      bz     c 


23      02         PC       23       P 
Using  F,   =  F,     +  F^  ^^  and  A^     =  C^     +  A^   , 

^  bz     bz      bz         bz       bz       bz  ' 


we  can  write: 

(V-  1)F,   +L  -  A^  ^^ 

b    =  ^ ^ ^ (2  12) 

bz         23       23       *  v^.x^; 

^        (A     +  A,    ■^) 
c       bz   ' 

Since  the  thermal  power  of  an  operating  reactor  is  essen- 

23      23 

tially  constant,  then  A    + A,    ,  the  total  rate  of  fissile 
^  c      bz 

consumption  in  the  reactor,  can  be  treated  as  a  fixed  value. 
It  can  be  seen  then,  that  Eq .  (2.12)  expresses  the  dependence 
of  the  rate  of  U233  conversion  in  the  breed  zone  on: 

(a)  the  rate  of  neutron  leakage  from  the  burn  region 
into  the  breed  zone,  [L  ] , 

(b)  the  rate  of  neutron  production  in  the  breed  zone, 
[(v-  1)  F^^]  and 

PC 

(c)  the  rate  of  neutron  loss  in  the  breed  zone,  [A,    ] • 

Losses  due  to  parasitic  absorption  and  leakage  in  the 
breed  zone  have  the  effect  of  cancelling  the  gains  made  by 
the  production  of  neutrons  in  the  zone.   It  is  these  factors 
that  must  be  reduced  while  simultaneously  increasing  the  rate 
of  neutron  production  in  the  breed  zone  to  secure  high  neutron 
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availability  and  improve  the  rate  of  U233  conversion. 
Combining  Eqs.  (2.10)  and  (2.11)  yields: 

C  °2       02 
b^  +  b,   =  ^c      ^bz 
""     ^^    ~~23 ^T  (2.13) 

c       bz 

which,  of  course,  is  the  system  breeding  ratio,  expressed  in 
I  Eq.  (2.1).   It  is  clear  then,  that  the  rate  of  breeding  in  the 
breed  zone  is  increased  only  at  the  expense  of  breeding  in  the 
burn  zone.   The  objective  then  becomes  to  optimize  the  total 
'  rate  of  U233  conversion  for  the  system,  under  the  constraints 
that  criticality  be  maintained  —  which  will  be  made  more  diffi- 
cult if  the  burn  zone  breeding  ratio  is  small.   Since  this  breed/ 
j  burn  concept  is  designed  for  application  in  commercial  power 
plants,  high  burnups  and  good  ore  utilization  are  essential. 
Design  parameters  that  will  allow  the  production  of  U233  in  the 
required  quantities  while  achieving  reasonable  fuel  element 
lifetimes  and  core  power  densities  must,  therefore  be  considered. 

2 . 5   Evaluation  of  Homogeneous  Zone  Characteristics 

Operation  of  the  breed/burn  fuel  cycle  for  the  pebble  bed 
reactor  system  has  been  defined  and  its  operation  described  in  the 
preceding  sections.   To  predict  more  precisely  the  characteristics 
of  such  a  system,  the  parameters  involved  in  breed/burn  operation 
for  homogeneous  one-zone  regions  will  be  evaluated.   This  analysis 
should  show,  for  a  given  composition,  how  a  region  would  operate 
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without  the  coupling  effects  (leakage)  of  a  zone  of  different 
composition  on  its  boundaries. 

Neutron  economy  in  any  reactor  system  is  a  function  of 
the  ambient  neutron  spectrum,  axial  and  radial  zoning  of 
compositions,  power  level  and  system  geometry.   The  rate  of 
neutron  leakage,  being  dependent  on  the  size  of  the  reactor, 
decreases  v/ith  increasing  core  volume,  holding  other  para- 
meters constant.   For  this  analysis,  to  provide  the  most 
ideal  and  hence  a  limiting-case  representation  of  the  feasi- 
bility of  breed/burn  operation,  an  infinite-medium  model  will 

be  used.   Pebble  bed  reactors,  designed  in  cylindrical  geo- 

3 

metrics,  operate  at  low  power  densities  (typically  5  w/cm 

[W3] )  and  hence  are  inherently  large  systems.   Reactors  having 

thermal  power  outputs  currently  competitive  commercially 

(^3000  MW,,  )  would  therefore  have  core  volumes  on  the  order  of 
th 
3 
600  m  ,  and  in  such  cases  would  have  neutron  leakage  fractions 

below  5%  [W3] .   Using  an  infinite  system  model  for  this  analy- 
sis, therefore,  is  not  an  unreasonable  approximation. 

For  layered  zoning  of  the  pebble  bed  system  described  in 
Chapter  1,  both  breed  and  burn  zones  will  be  treated  here  as 
uniform,  separate  regions.   In  this  way,  for  a  given  spectrum, 
the  best  possible  operating  characteristics  can  be  determined 
as  there  will  be  no  coupling  effects  between  regions  which 
would  perturb  the  spectrum  in  the  vicinity  of  the  boundary. 
As  shov7n  in  the  previous  section,  the  breeding  ratio  for  the 
system  is  simply  the  sum  of  the  breeding  ratios  in  the  breed 
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and  burn  zones.   Therefore,  if  one  can  determine  the  maximum 
amount  of  breeding  possible  in  a  homogenized  region,  this  will 
approximate  the  system-averaged  breeding  ratio  for  a  two-zone 
reactor  in  the  ideal  case. 

The  pebble  bed  HTGR  utilizes  injection-molded  spherical 
graphite  elements,  illustrated  in  Fig.  [2.1],  fueled  with 
pyrolytic  carbon  (PyC)  and  silicon  carbide  (SiC)  coated  ThO^ 
and  UO2  microspheres.   The  balls  themselves  provide  all  of 
the  fuel,  moderator,  structural  material  and  coolant  channels 
(free  space  between  balls)  necessary  for  the  reactor.   There- 
fore, in  a  clean  core,  the  only  materials  present  are  those 
found  in  the  newly  fabricated  feed  balls,  in  addition  to  the 
helium  coolant.   As  a  function  of  the  carbon- to-heavy  metal 
atom  ratio,  the  number  densities  of  fertile  feed  ball  materials 
are  given  in  Table  [2.1]  (for  calculation  of  number  densities, 
refer  to  Appendix  D) . 

Variation  in  the  carbon-to-heavy  metal  atom  ratio  signifi- 
cantly alters  the  neutron  spectrum  in  a  given  region  due  to  the 
change  in  moderation.   One-group  constants  for  carbon-to-heavy 
metal  atom  ratios  of  110,  250,  325  and  450  are  given  in  Appen- 
dix A,  collapsed  from  the  AMPX-generated  4-group  cross-sections 
given  in  Appendix  C  (the  AMPX  system  is  described  in  Chapter  3) , 
To  illustrate  the  effect  of  the  carbon-to-heavy  metal  atom 

ratio  on  the  neutron  spectrum,  the  spectrum-averaged  values 

23  •  • 

for  n   ,  the  rate  of  neutron  production  per  fissile  absorption, 

given  by: 
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FIGUP^  2.1:   PEBBLE  BED  HTGR  FUEL  ELEMENT 
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TABLE  2.1 


COMPOSITION*  OF  PEBBLE  BED  FEED  BALLS 


C/HM  ATOM  RATIO 
NUCLIDE       110        250        325        450 

Thorium    2.422E-4   1.041E-4   7.872E-5  5.702E-5 

Oxygen     4.844E-4   2.083E-4   1.574E-4  1.140E-4 

Carbon     2.666E-2   2.610E-2   2.600E-2  2.591E-2 

Silicon    6.644E-4   2.857E-4   2.160E-4  1.564E-4 


*  Number  densities  in  b   cm 
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23 

23    ^°f 

n   =  oT  /  (2.14) 

a   -^ 

a 

the  microscopic  fission  cross-section  for  Th232,  o^^^,  and  the 

microscopic  neutron  absorption  cross-sections  for  the  parasitic 

w    V.        carbon    oxygen    .,    silicon        ,  ,    .  . 
absorbers,  a^       ,  o^      ^^        and  o^        j-^v^xi^  ^^^  tabulated  in 

Table  [2.2] . 

It  can  be  seen  from  these  results  that  the  preference  for  a 

hardened  neutron  spectrum  is  not  immediately  justified.   The 

values  given  in  Table  [2.2]  are  those  used  in  determining  the 

breeding  ratio  of  the  system,  expressed  by  Eq.  (2.1).   It  is 

immediately  apparent  that  increasing  the  carbon-to-heavy  metal 

23 
atom  ratio  from  110  to  4  50  increases  n   .   This  is  opposed  to 

23 

our  analysis  m  that  an  increase  in  n    with  added  moderation 

would  cause  an  increase  in  the  breeding  ratio.   However,  Table 
[2.2]  also  reveals  decreases  with  a  harder  spectrum  in  the 
absorption  cross-sections  of  the  parasitic  nuclides  and  an 
increase  in  the  fission  cross-section  of  Th232.   These  factors 

lead  to  increases  in  the  breeding  ratio  v/ith  spectrum  harden- 

23 
ing.   In  fact,  while  n    decreases  ^2%  from  C/HM  =  450  to 

C/HM  =  110,  the  macroscopic  parasitic  absorption  cross-section 

(Z  ^  +  Z  °  +  E  ^"^)  decreases  ^^^36%,  v/hile  the  Th232  fission 
a      a      a 

cross-section  increases  'v72%.   So  although  Eq.  (2.1)  is  domi- 

23  •   • 

nated  by  n   ,  the  much  larger  changes  in  Th232  fissioning 

and  parasitic  absorption  could  have  the  dominant  effect  on 

the  breeding  ratio.   For  these  reasons  it  is  important  to 

more  carefully  consider  the  effects  of  the  ambient  spectrum 

on  breed/burn  operation. 
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TABLE  2 . 2 
SPECTRUM-AVERAGED  HOMOGENEOUS  ZONE  PARAMETERS 

C/HM  ATOM  RATIO 
PARAMETER       110        250        325        450 
n^^  2.209      2.240      2.246      2.253 


a^2        8.097E-3   6.051E-3   5.446E-3   4.718E-3 

a^  3.830E-4   7.095E-4   8.191E-4   9.507E-4 

a 

o^  5.306E-4   4.146E-4   3.811E-4   3.409E-4 

a 

O^^  1.358E-2   2.709E-2   3.160E-2   3.702E-2 

a 


Cross-sections  in  barns 
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2 . 6   Variation  in  Neutron  Spectra 

The  proposal  by  Fischer  et  al  [Fl]  for  implementation  of 
the  breed/burn  fuel  cycle  in  a  hard  spectrum  core  indicated 
that,  with  respect  to  breeding  ratio,  the  best  results  would 
be  obtained  if  the  ambient  spectrum  was  hardened  as  much  as 
possible.   In  considering  the  pebble  bed  HTGR,  our  first  indi- 
cations verify  Fischer's  claims  for  the  preference  of  the 
hardened  spectrum.   It  has  now  become  necessary  to  justify 
this  claim  more  fully  if  it  is  to  be  used  as  a  design  objec- 
tive for  thermal  reactor  versions  of  the  breed/burn  concept. 

Two  criteria  for  successful  operation  of  the  breed/burn 
cycle  have  been  established:   (1)  the  ability  to  maintain 
system  criticality,  and  (2)  the  ability  to  breed  adequate 
fissile  mass  in  the  breed  zone  to  fuel  the  burn  zone.   The 
second  objective  is  solely  a  function  of  the  breeding  ratio  in 
the  breed  zone.   Maximizing  the  rate  of  U233  conversion  by 
increasing  the  parameters  in  Eq .  (2.1),  as  described  in  the 
previous  section,  makes  it  easier  to  achieve  the  second  goal 
alone,  if  adequate  exposure  time  can  be  provided  by  artifi- 
cially satisfying  the  first  objective.   However,  when  the 
criteria  for  successful  operation  must  be  satisfied  simul- 
taneously, it  is  realized  that  the  optimization  of  one 
parameter  detracts  from  the  ability  to  satisfy  the  other. 

By  this  reasoning,  Fischer  et  al,  in  achieving  both 
breed/burn  operation  and  the  highest  breeding  ratio  possible, 
might  not  be  able  to  maintain  criticality  in  their  overall 
system.   If  this  is  the  case,  operating  a  breed/burn  cycle 


-42- 


becomes  a  futile  task.  Then,  not  only  must  the  factors  in- 
volved in  obtaining  a  high  breeding  ratio  be  considered,  but 
also  those  factors  affecting  the  criticality  of  the  system. 
Sheaffer  et  al  [SI]  have  correlated  microscopic  cross- 
sections  with  the  ambient  neutron  spectrum  in  fast  reactors 
by  defining  a  spectral  index  param.eter,  S,  given  by: 


(2.15) 


f      tr 


Physically,  S  is  a  measure  of  the  average  neutron  energy  in 
the  core  population.   Cross-sections,  as  a  function  of  S,  are 
given  by: 


a  =  a,S^  (2.16) 


where  o -.    and  g  are  constants  appropriate  for  the  indicated 
nuclide  and  cross-section  (the  constants  a,  and  g  are  tabula- 
ted in  Appendix  B,  Table  [B.l]  for  the  materials  used  in  this 
analysis).   The  correlation  expressed  by  Eq .  (2.16)  was  devel- 
oped for  application  to  fast  reactor  systems  and  has  been 
verified  (v;ith  very  good  accuracy)  for  values  of  S  in  the 
range  0.3  to  1.0  [SI]  (corresponding  to  a  wide  range  of  fast 
reactor  spectra,  from  oxide-fueled  LMFBRs  to  metal-fueled 
GCFRs)  .   Sheaffer 's  method  can  be  used  to  evaluate  breed/burn 
performance  in  FBR  systems. 
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The  ability  of  the  system  to  breed  fissile  enrichments 
suitable  to  fuel  its  fast  core  must  be  investigated.   In  any 
breeding  system,  if  it  was  assumed  that  the  neutron  flux  could 
be  maintained  artificially,  then  the  system  would  eventually 
achieve  a  conversion  ratio  of  1.00,  at  which  time  the  rate 
of  fertile  neutron  capture  would  exactly  equal  the  rate  of 
fissile  destruction  by  neutron  absorption.   Thus,  the  system 
would  maintain  an  asymptotic  fissile  enrichment  limit  which 
would  be  the  maximum  possible  enrichment  compatible  with  the 
am>bient  spectrum. 

If  the  fertile  capture  rate  is  equal  to  the  fissile  ab- 
sorption rate: 


No   (f) 
c 


a^ 


fertile 


(2.17) 


fissile 


Neglecting  all  heavy  metals  except  the  principal  fertile  and 
fissile  species,  the  fissile  enrichment,  e,  can  be  defined  as: 


N 


fissile 


,  fissile  .  ^^fertile 
N         +  N 


(2.18) 


Combining  Eqs .  (2.17)  and  (2.18),  an  expression  for  the 

asymptotic  fissile  enrichment,  ^n^y'    ^^    ^    function  of  the 

, .     fertile,   fissile  . 
cross-section  ratio  a         /a  is: 

c         a 


'ASY 


fertile  ,   fissile 
a        /a 
c ^  a 

/   fertile,   fissile 
1  +  (a        /a^ 


)   ■ 


(2.19) 
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For   a   plutonium-239/uraiiium-238    system,    Eq .     (2.19)    becomes: 

a        /a 

^ASY   "    ,    ^     .       fr;      T97  (2.20) 

1   +    (o        /a        ) 
c      ^     a 


and  similarly  for  a  U233/Th  system; 


02  ,   23 

a    /a 

e     =  2 ^ (2  21) 

ASY    .  ,  ,   02/   23,   .  ^^.^J-J 

1  +  (a    /a    ) 
c     z 


With  these  approximations  the  practical  end-of-life  (EOL) 
fissile  enrichment,  ^^qt /  obtainable  with  high  fuel  burnup 
in  a  breed/burn  system,  can  be  roughly  estimated. 

Consider  the  following  conceptual  picture  of  breed/burn 
operation:   pairs  of  physically  identical  fuel  elements  exist 
—  one  in  the  breed  zone  and  the  other  in  the  burn  zone.   The 
fissile  enrichment  in  the  breed  zone  at  end-of-cycle  must 
equal  the  fissile  enrichment  in  the  burn  zone  at  the  beginning- 
of-cycle  if  a  "stand-alone"  breed/burn  cycle  is  to  be  feasible. 
For  every  fission  in  the  burn  zone  element  roughly  one  excess 
neutron  is  produced  which  can  be  used  to  produce  fissile 
material  in  the  breed  zone  element.   Since  both  elements  have 
very  nearly  equal  heavy  metal  contents,  the  maximum  percent 
fissile  enrichment  in  the  breed  element  is  equal  to  the  percent 
burnup  in  the  burn  zone.   That  is,  for  breed/burn  operation 


(in  %)  : 


B  [MWD/MTHM] 

^BOL^  =  ^EOK    =^      ,^4 ^2.22) 

C        bz        10       • 
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The  potential  of  different  reactor  types  is  as  follows: 

LWR:   ^bOl'^-^^'  ^EOl'^^^'^^^    MWD/MTHM=  3% 

HTGR:   £gQ^%9%;  Bg^^%120,000   MWD/MTHM=  12% 

LMFBR  or  GCFR:  £^^^'^15%;  B^^^ 'V  100 ,  000   MWD/MTHM  =  10%. 

This  information  is  illustrated  in  Fig.  [2.2],  v/here  it  can 
be  seen  that  the  potentially  feasible  breed/burn  systems  would 
be  those  operating  below  the  line,  whereas  those  above  the  line 
would  not  be  feasible.   The  circumstances  indicated  in  Fig. 
[2.2]  are,  in  fact,  a  large  part  of  the  motivation  for  investi- 
gating the  HTGR  as  a  breed/burn  reactor  in  the  present  study. 

Since  the  cross-sections  of  the  fertile  and  fissile 
species  vary  with  the  ambient  spectrum,  Eq .  (2.19)  can  be 
used  to  estimate  the  attainable  asymptotic  enrichment  limit 
for  any  given  system.   This  limit  indicates  for  a  reactor 
operating  with  the  breed/burn  fuel  cycle,  whether  to  harden  or 
soften  the  ambient  spectrum.   A  high  asymptotic  enrichment 
would   imply  that  relatively  short  burnups  would  be  required 
to  achieve  a  desired  enrichment  value,  whereas  a  low  asymptotic 
limit  would  mean  that  if  a  beginning-of -lif e  enrichment  was 
required  near  the  asymptotic  value,  then  very  high  burnups 
would  be  essential  to  achieve  that  limit.   This  argument  adds 
preference  to  those  spectra  whose  asymptotic  enrichment  limits 
are  as  high  as  possible. 

From  Eq.  (2.16),  cross-sections  can  be  deterroined  for  any 
given  fast  reactor  spectrum  (composition) .   For  this  analysis. 
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FIGURE  2.2:   SCHEMATIC  DIAGRAM  SHOWING  POTENTIAL 
FEASIBILITY  OF  BREED/BURN  SYSTEMS 


•Typical  FBR  with  current  burnup  limits. 
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the  following  approximation  is  used 


^  fertile  ^  ^  fertile  (2.23) 

c  a 


which,  from  Eq.  (2.19),  yields: 

fertile,   fissile 

^a       /^a 

^ASY  ""   7~~~,      fertile,   fissile,  *  (2.24) 

1  +  (o__       /a        ) 
a       'a 

This  approximation  overestimates  the  value  of  e  and  there- 

fore  yields  a  conservative  upper  limit  for  the  asymptotic 
enrichment.   Listed  in  Table  [2.3]  are  values  for  £,c:v'  using 
Eq .  (2.24),  calculated  for  spectra  varying  from  S=0.3  to 
S=1.0.   Also  tabulated  are  values  for  fissile  n ,  given  by: 


n  =  -^—  (2.25) 

a 


for  the  indicated  fissile  species.   It  is  clear  from  the  data 
in  Table  [2.3]  that  values  for  fissile  n  for  both  the  Pu239/U 
and  U233/Th  systems  increase  with  spectrum  hardening,  justify- 
ing the  claim  of  Fischer  et  al  that  harder  spectra  will  yield 
better  breeding  potential.   This  alone,  however,  does  not 
justify  the  preference  of  the  harder  spectra  for  the  breed/ 
burn  concept.   In  fact,  because  a  balance  must  exist  between 
fissile  breeding  and  system  criticality,  the  asymptotic  fissile 
enrichments,  illustrated  in  Fig.  [2.3],  indicate  that  the 
softer  fast  reactor  spectra  would  be  preferable.   Operation 
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TABLE  2 . 3 

ASYMPTOTIC  ENRICHMENTS  AND  FISSILE  n  VALUES 

FOR  FAST  SPECTRA* 


23  23  ^9  k^ 

ASY      "  ASY       ' 


0.3  0.1388  2.175  0.1480  2.308 

0.4  0.1210  2.224  0.1368  2.415 

0.5  0.1085  2.264  0.1286  2.502 

0.6  0.09918  2.296  0.1223  2.574 

0.7  0.09190  2.324  0.1172  2.638 

0.8  0.08591  2.348  0.1128  2.694 

0.9  0.08097  2.370  0.1091  2.744 

1.0  0.07675  2.390  0.1059  2.790 


*  From  Sheaf fer  et  al  [SI]  correlation:   cross-sections 
are  listed  in  Table  [B.2]. 
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j   can,  in  this  case,  achieve  much  higher  fissile  enrichments  for 
a  given  burnup  than  can  a  system  operating  in  a  much  harder 
spectrum.   For  a  Pu239/U  fueled  system,  the  softer  (S=0.3) 
spectrum  burnup  (assuming  a  linear  relation  between  burn  zone 
and  breed  zone  enrichments)  would  require  an  exposure  time  ^ 
32%  less  than  the  hard  (S=1.0)  spectrum,  while  for  the  U233/ 
Th  system  the  same  change  in  spectrum  would  reduce  the  exposure 
time  by  '^^  4  8%. 

These  cases  have  indicated,  contrary  to  the  proposal  of 
Fischer  et  al,  that  a  softer  spectrum  is  preferred  for  breed/ 
burn  operation  in  a  fast  reactor.   A  similar  analysis  can  be 
performed  to  determine  the  more  favorable  operating  conditions 
for  the  thermal/epithermal  HTGR.   Because  Sheaffer's  correla- 
tions were  developed  for  fast  reactor  calculations  they  do  not 
account  for  the  high  degree  of  resonance  shielding  in  the 
fertile  species  for  highly  epitherraal  applications.   Hence, 
the  model  is  invalid  for  values  of  S<  0.3;  the  lowest  values 
verified. 

A  similar  spectral  index  can  be  derived  for  application 
in  the  U233/Th  HTGR  cases.   In  the  expression  for  the  reso- 
nance escape  probability  from  classical  continuous  slowing 
theory,  p,  given  by: 


exp 


E   dE' 
a 


J 


E   E'  U^^  +  Z3 


(2.26) 
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an  effective  Z^  can  be  used  to  account  for  leakage  by  adding 
a  "leakage  cross-section": 


Z  ->  E   +  DB  . 
a    a 


(2.27) 


Then,  accounting  for  criticality  by  applying  the  neutron 
balance: 


Z^  +  DB   =  vZ^   , 
a  r 


(2.28) 


one  gets: 


p  -  exp 


; 

E 


vZ 


^(vZ^  +  Eg) 


dE' 
E' 


(2.29) 


Since  Z   -  Z    ,  and  recognizing  that  in  the  thermal  /  epither- 
mal  region  vZ^ .  «  Z    ,  a  spectral  index,  S'  can  be  defined 
by: 


vZ 


S' 


^^ 


(2.30) 


tr 


Using  this  definition,  the  index  can  be  redefined  as: 


S' 


1  +  S 


I   / 


(2.31) 


which  yields  the  same  form  used  by  Sheaffer  et  al  [SI] ,  given 

by: 

vZ^ 

(2.15) 


vZ^  +  C2. 
f    ^  tr 
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From  the  computations  to  be  discussed  in  Chapter  3,  values 
used  to  calculate  S  were  found  for  the  pebble  bed  HTGR  composi- 
tions used  in  the  analysis  of  Section  [2.5].   Table  [2.4]  lists 
values  for  ^^gy'  9^^®^  ^Y  Eq .  (2.21),  with  corresponding  val- 
ues of  S  as  a  function  of  the  carbon- to-heavy  metal  atom  ratio. 

2  3 
In  Fig.  [2.4]  the  asymptotic  enrichments,  e,ev   '  ^^^  plotted 

versus  the  spectral  index  for  both  the  U233/Th  fast  reactor 
cases  and  the  pebble  bed  systems.   It  can  be  seen  that  in  the 
thermal/epithermal  region,  spectral  hardening  is  advantageous 
in  increasing  the  asymptotic  enrichment  limit. 

As  shown  previously,  although  the  harder  HTGR  spectra 
lead  to  increasing  values  of  fissile  n^  decreasing  parasitic 
absorption  rates  and  improved  fertile  fission  rates  combined 
with  the  benefit  of  higher  asymptotic  enrichment  limits  indi- 
cate a  strong  preference  for  the  use  of  a  hardened  spectrum 
for  successful  thermal/epithermal  breed/burn  operation.   This 
conclusion  is  contrary  to  that  of  fast  reactors,  where  the 
softer  spectra  are  preferred.   The  indication  here  is  that 
an  optimum  breed  zone  ambient  spectrum  may  be  attainable  in 
the  regions  between  those  considered  in  this  analysis.   The 
apparent  maximum  value  for  e         ,    shown  in  Fig.  [2.4],  clearly 
indicates  this  possibility.   Hence,  in  order  to  maximize  the 
effectiveness  of  the  breed/burn  concept  it  would  be  advisable 
to  consider  systems  operating  with  neutron  spectra  in  the 
area  of  maximum  ^j,qy   between  a  "soft"  FBR  and  a  "hard"  HTGR. 
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TABLE    2 . 4 
SPECTRAL    INDEX    AND    ASYMPTOTIC    ENRICHMENT 
LIMITS    FOR   THE    U233/Th   PEBBLE    BED    HTGR 

C/HM   ATOM    RATIO 
PARAMETER  110  250  325  450 

S  0.0724         0.0583         0.0518         0.0431 


e  0.03526      0.02686      0.02497      0.02302 

/\0  jl 
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2 . 7   System  Criticality 

In  the  previous  section  it  was  shown  that  in  using  a 
breed/burn  fuel  cycle  in  the  pebble  bed  HTGR,  it  is  advanta- 
geous to  harden  the  breed  zone  neutron  spectrum  in  order  to 
increase  the  rate  of  fissile  U233  production.   The  required 
burn  zone  burnup  for  each  element  necessary  to  achieve  any 
given  breed  zone  fissile  enrichment  can  therefore  be  shorten- 
ed, and  hence,  the  number  of  passes  through  the  breed  zone 
required  of  each  fuel  ball  can  be  minimized. 

The  fissile  enrichments  attainable  for  each  composition 
can  be  approximated  using  Eq.  (2.21),  again  utilizing  the 
one-group  constants  given  in  Appendix  A.   As  a  function  of 
the  carbon- to-heavy  metal  atom  ratio,  these  enrichments  are 

listed  in  Table  [2.5]  along  with  the  cross-section  ratio, 

02    23 
a    /a    .   Using  these  values  to  approximate  the  beginning- 
c     a 

of-life  fissile  enrichment  for  a  clean  core,  the  system  criti- 
cality at  the  beginning  of  each  cycle  can  be  estimated.   Apply- 
ing a  one-group  model,  the  infinite-medium  multiplication 
factor,  k   ,  can  be  estimated  at  the  beginning-of-lif e  as: 


k      -  ^  (2.31) 

BOL     a 


and  the  BOL  excess  reactivity,  P^qt  '  ^^ • 

k      -  1 
"boL  (2.32) 

^BOL      k 

CO 

BOL 
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TABLE  2.5 
ASYMPTOTIC  FISSILE  ENRICHMENT  LIMITS 

C/HM  ATOM  RATIO 
PARAMETER    110       250       325       450 

'23      0.03655   0.02760   0.02561   0.02357 
a 

e  (%)■         3.526  2.686         2.497  2.302 
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Values  for  k^     and  p     are  given  in  Table  [2.6]. 
BOL       ^^^ 
It  can  be  seen  that  in  all  cases,  k^     >  1.00,  indicating 

°°BOL 
criticality  at  the  start  of  the  cycle.   Using  the  amount  of 

U233  bred,  for  example,  in  the  region  of  C/HM  =  110  to  fuel 

the  regions  where  C/HM  =  250  or  higher  would  yield  an  even 

greater  value  for  k^     .   For  the  four  compositions  consid- 

BOL 
ered,  the  values  of  k^     obtained  for  various  combinations  of 

BOL 
breed  and  burn  zone  carbon- to-heavy  metal  atom  ratios  are 

shown  in  Table  [2.7].   The  highest  k     obtained  is  for  that 

-^  CO 

BOL 
of  a  burn  zone  with  C/HM  =  450  being  fueled  by  a  breed  zone 

where  C/HM  =  110.   This  combination  has  the  softest  burn  zone 
spectrum  operating  with  the  hardest  spectrum  (and  hence,  high- 
est ^T^ay)    breed  zone  of  those  compositions  evaluated. 

Careful  analysis  is  required  in  these  cases  to  insure  that 
a  given  breed  zone  would  produce  not  only  adequate  fissile 
enrichment,  but  also  sufficient  total  fissile  mass  to  fuel  the 
burn  zone.   From  this  analysis,  the  first  order  estimate  leads 
to  the  conclusion  that  it  is  possible  to  breed  enough  fissile 
material  in  a  breed  zone,  given  adequate  exposure  time,  to  fuel 
and  obtain  criticality  in  a  pebble  bed  breed/burn  system.   The 
amount  of  BOL  excess  reactivity  is  quite  small,  however,  which 
indicates  that  it  may  not  be  possible  to  maintain  criticality 
long  enough  to  breed  fissile  enrichments  near  the  asymptotic 
limits.   One  can  adapt  to  this  situation  to  some  extent  by 
using  rapid  fuel  flow  through  the  core,  approaching  continuous 
refueling  in  the  limiting  case,  but  of  course  this  increases 
the  complexity  of  operation. 
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TABLE  2.6 
BEGINNING-OF-LIFE  INFINITE-MEDIUM 
MULTIPLICATION  FACTORS 

C/HM  ATOM  RATIO  (both  zones) 
PARAMETER     110        250        325        450 


vE^      1.435E-3   7.782E-4   6.163E-4   4.662E-4 


E        1.293E-3   7.110E-4   5.695E-4   4.392E-4 
a 


k   ^     1.110      1.095      1.082      1.062 
o°BOL 


p^^^  (%)   9.91       8.64       7.60       5.71 
BOL 
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TABLE  2 . 7 


k^  FOR  VARIOUS  ZONE  COMBINATIONS* 


BURN  ZONE  C/HM  ATOM  RATIO 
110      250      325      450 

110   1.110    1.256    1.286    1.313 
BREED  ZONE 

250   0.9515   1.095    1.125    1.152 
C/HM  ATOM 

325   0.9103   1.052    1.082    1.109 
RATIO 

450   0.8654   1.005    1.035    1.062 

*  k^  in  burn  zone  using  asymptotic  enrichments 
achievable  in  breed  zone 
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The  analysis  of  Section  [2.6]  can  be  refined  further  to 
consider  some  of  the  specific  features  of  the  breed/burn  cycle 
In  the  expression  for  k^  ,  given  by  Eq .  (2.31),  a  "leakage 
cross-section"  can  be  added  to  the  absorption  cross-section, 
so  that'. 


k  =  y — -^  (2.33) 

^a  ^  ^L 


2 
v;here,  for  example,  one  can  use  E   =  DB   .   The  criterion  for 

criticality,  that  k>^1.00,  yields  the  necessity  that: 

k^  1  1  +  ^  (2.34) 

a 
Z  23 

^oo  -  11-23— f-Or—rP  (2.35) 

a      a      a 

p 
where  E    refers  to  neutron  absorption  by  the  parasitic  ab- 
sorbers . 

Therefore,  Eq .  (2.35)  can  be  used  to  find  that: 


k     1^  Q2    23        P    23     (2.36) 


where  b   is  the  breeding  ratio  in  the  burn  zone,  approximated 
by: 

a 
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02,^  23_,  ,  P    23 
a  ''^  'a 

criticality  limit  becomes: 


In  Eq.  (2.37),  both  Y. ^"^ "■  / Z ^"^ " and    Z^V^.    are  small.   The 


^co  >  1  +  1  +  b    •  (2.38) 

c 


In  the  burn  zone,  however,  the  maximum  value  of  b   is  attained 

c 

under  EOL  asymptotic  conditions.   In  this  case,  b   =1.0, 
which  therefore  defines  a  lower  limit  for  k   of: 


k^  t   1-5  .  (2.39) 

which  is  now  a  third  criteria  for  successful  system  operation. 

Thus,  k^  for  the  burn  zone  must  exceed  1.5  for  a  breed/burn 

cycle  to  be  able  to  operate  independent  of  additional  external 

fissile  sources. 

Since  the  greatest  value  obtained  for  k      ,  given  in 

BOL 
Table  [2.7],  was  k^  "^   1.3,  we  must  conclude  that  for  the  pebble 

bed  HTGR,  stand-alone  breed/burn  operation  does  not  appear 

possible. 

2 . 8   Summary 

In  this  chapter  the  parameters  affecting  successful  breed/ 
burn  operation  in  the  pebble  bed  HTGR  have  been  examined.   While 
preliminary  evaluation  suggested  that  spectral  hardening  would 
detract  from  system  performance  due  to  the  decreasing  values  of 
fissile  T\ ,    closer  examination  revealed  that  while  n  decreases 
slightly,  parasitic  absorption  decreases  significantly  while 
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the  fertile  fission  rate  also  increases  substantially.   In 
addition,  harder  spectra  yield  higher  asymptotic  enrichment 
values,  thus  making  it  easier  to  achieve  acceptable  fissile 
enrichments  in  the  breed  zone  with  reasonable  burnups  in  the 
burn  zone.   The  results  also  shov/ed  that  highly  thermalized 
systems  are  best  for  the  burn  zone  because  of  their  high  k^ 
at  a  given  load  enrichment  and  low  asymptotic  (hence  discharge) 
fissile  enrichment.   Of  all  reactors  of  current  interest, 
it  v/as  found  that  the  HTGR  had  the  most  favorable  ratio  of 
discharge  burnup  ro  reload  enrichment. 

These  observations,  with  respect  to  spectral  hardening, 
were  found  to  apply  for  thermal  and  epithermal  systems  only. 
In  fact,  for  fast  reactors  it  was  found  that  the  breed  function 
of  breed/burn  cycles  would  be  best  carried  out  in  softer 
spectra,  such  as  those  found  in  oxide-fueled  LMFBRs  --  or, 
indeed,  even  softer. 

In  examining  system  criticality,  although  it  was  found 
that  a  breed  zone  in  a  pebble  bed  reactor  could  fuel  and 
achieve  criticality  (k   =1.0)  in  an  infinite  burn  zone,  the 
high  BOL  excess  reactivity  (k^  >_   1.5)  necessary  for  finite 
coupled  zone  (high  burn  and  breed  zone  leakage)  system  opera- 
tion could  apparently  not  be  achieved.   Hence,  operating  such 
a  system  in  the  pebble  bed  HTGR  would  require  some  degree  of 
fissile  makeup  to  operate  for  any  sustained  period;  either 
that  or  a  substantially  harder  breed  zone  spectrum  (yielding 
higher  asymptotic  enrichment  and  hence  larger  k^  in  the  burn 
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zone) .   An  even  higher  carbon-to-heavy  metal  atom  ratio  in  the 
burn  zone  could  also  be  contemplated. 
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III.   COMPUTATIONAL  METHODS 

3 . 1  Introduction 

The  analysis  in  the  previous  chapter  led  to  the  conclusion 
that  in  a  pebble  bed  HTGR,  it  would  probably  not  be  possible 
to  sustain  a  stand-alone  breed/burn  system.   However,  the  one- 
group  calculations  performed  could  only  indicate  the  gross 
features  of  system  operation  and  were  limited  to  either  the 
beginning-of-lif e  or  the  asymptotic  limit  of  ultra-high  burnup. 
The  effects  of  fission  product  buildup,  U235  and  Pa233  produc- 
tion, Th232  depletion,  fast  fission,  power  density  and  group- 
wise  changes  in  flux  and  cross-sections,  all  as  functions  of 
burnup,  were  neglected.   When  fully  evaluated,  the  sum  of  the 
above  effects  can  significantly  alter  the  reactor  performance 
characteristics  as  a  function  of  time,  which  therefore  neces- 
sitates a  more  precise  evaluation. 

Calculation  of  the  performance  characteristics  of  the 
proposed  breed/burn  system,  using  state-of-the-art  computer 
programs,  is  described  in  this  chapter.   The  system  model, 
cross-section  generation,  nuclear  data  input,  codes  used  and 
problem  descriptions  are  described. 

3 . 2  System  Model 

In  the  breed/burn  fuel  cycle,  it  is  necessary  for  the 
reactor  system  to  operate  under  "breakeven"  breeding  condi- 
tions.  That  is,  as  a  minimum,  one  fissile  atom  must  be  pro- 
duced in  the  breed  zone  for  every  one  destroyed  in  the  burn 
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zone.   A  first  approximation  can  be  used  to  evaluate  this 
ability  (i.e.,  by  considering  a  homogeneous  system)  as  v;as 
done  in  Chapter  2.   The  two-zone  layered  system,  described  in 
Chapter  1,  will  operate  with  a  strong  neutronic  coupling  be- 
tween zones.   The  burn  zone  fission  rate,  and  hence,  power 
density,  will  be  higher  than  that  in  the  breed  zone.   Also, 
the  total  neutron  flux  in  the  burn  zone  should  be  higher, 
whereas  groupwise  flux  will  be  shifted  towards  the  faster  groups 
in  the  breed  zone,  assuming  that  the  carbon-to-heavy  metal  atom 
ratio   in  the  breed  zone  is  less  than  that  of  the  burn  zone. 

A  conservative  overestimation  of  the  performance  of  the 
system  can  be  found,  therefore,  by  homogenizing  the  two  zones 
into  a  single  region.   If  the  conditions  required  for  success- 
ful breed/burn  operation  can  be  met  in  the  homogeneous  region, 
then  a  more  careful  analysis  is  required  to  determine,  when  the 
reactor  is  modeled  in  a  heterogeneous,  axially-layered  configura- 
tion, if  satisfactory  performance  is  still  achieved. 

In  analyzing  the  heterogeneous  two-zone  system,  the  re- 
quirement for  the  breeding  of  adequate  total  mass  in  the  breed 
zone  to  fuel  the  burn  zone  dictates  size  limitations  on  the  two 
zones,  contingent  upon  the  number  of  passes  through  each  zone 
required  to  achieve  design  burnup  limits.   That  is,  given  the 
compositions  of  each  zone,  in  a  period  of  one  cycle,  the  total 
fissile  mass  being  bred  in  the  breed  and  burn  zones  must  equal 
the  total  mass  destroyed  in  the  burn  zone  for  net  breeding; 
for  breed/burn  operation  the  breed  zone  alone  must  provide  the 
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burn  zone  fissile  deficit.  Therefore,  calculations  must  be 
performed  iterating  between  size  options  for  the  two  zones. 
If  the  reactor  is  assumed  to  be  infinite  radially,  then  the 
calculations  can  be  performed  using  an  infinite  slab  model 
varying  only  the  compositions  and  axial  thicknesses  of  each 
zone. 

The  first  analysis  to  be  performed,  however,  utilizes  a 
homogeneous  zone  as  previously  described.   Unity  albedo  boun- 
dary conditions  are  used  for  all  three  dimensions,  thus  effec- 
ting an  infinite  reactor  calculation.   Again,  an  infinite  slab 
treatment  is  appropriate.   As  will  be  described  in  Section  [3.4], 
compositions  (carbon-to-heavy  metal  atom  ratios  and  fissile 
enrichments)  will  be  varied  to  determine  system  characteristics. 

3. 3   Computer  Codes 

Analysis  of  the  pebble  bed  HTGR  breed/burn  fuel  cycle  in 
the  homogeneous  system  model  described  in  the  previous  section 
required  precise  definition  of  nuclear  data  parameters  for 
input  into  a  code  that  would  treat  the  one-dimensional  diffu- 
sion-depletion problem  accurately.   In  this  section  the  com- 
puter codes  used  for  this  application  and  their  respective 
input  parameters  are  described. 

3.3.1   PDQ-7 

The  PDQ-7  code  [Cl]  was  used  in  the  present  work  for  the 
calculation  of  the  operating  characteristics  of  the  breed/burn 
system  model.   PDQ-7  can  solve  the  diffusion  and  depletion 
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equations  in  up  to  three  dimensions  in  rectangular,  cylindri- 
cal, spherical  or  hexagonal  geometry  for  up  to  five  energy 
groups.   Depletion  calculations  and  management  of  the  neutron 
cross-section  data  are  handled  by  the  HARMONY  system  [B2], 
which  provides  a  flexible  representation  of  the  time-dependent 
depletion  equations  for  any  nuclide  chain. 

The  PDQ-7  code  solves  the  diffusion  equation  by  discre- 
tizing  the  energy  variable  and  finite  differencing  (central) 
the  equation  in  space.   For  this  application,  only  one-dimen- 
sional problems  were  treated,  which  the  code  solves  by  Gaussian 
elimination. 

For  all  calculations,  PDQ-7  requires  that  input  parameters 
define  the  system  geometry  and  energy  groups  with  respective 
chi  (fission  neutron)  spectra,  as  well  as  the  complete  repre- 
sentation of  nuclide  chains,  decay  constants,  initial  nuclide 
number  densities  and  microscopic  neutron  cross-section  data. 

The  problems  described  in  Section  [3.4]  were  all  solved 
using  four  energy  groups  in  an  infinite  slab  system,  with  de- 
pletion parameters  calculated  for  a  static  core  arrangement, 
from  the  beginning-of-life  to  an  end-of -cycle  fuel  burnup  of 
120, 000  MWD/MTHM. 

The  depletion  chain  treatment  followed  46  nuclides,  for 
which  the  respective  nuclide  chains  are  detailed  in  Appendix 
LD.5].   Initial  nuclide  number  densitites  were  calculated  as 
a  function  of  the  carbon- to-heavy  metal  atom  ratio  and  the 
beginning-of-life  fissile  enrichment,  as  described  in  Appendix 
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[D.3].   Beginning-of-life  heavy  metal  number  densities  are 
given  in  Table  [D.4]  as  a  function  of  fissile  enrichment,  and 
all  other  initial  concentrations  are  listed  in  Table  [D.5]. 

Exposure  times  (in  hours)  used  for  the  depletion  incre- 
ments were  calculated  for  burnup  intervals  as  detailed  in 
Appendix  [D.4].   All  problems  were  treated  for  burnups  to 
120, 000  MWD/MTHM. 

Cross-section  input  for  PDQ-7  was  provided  by  calculations 
using  the  AMPX  Code  System  [G2],  as  described  in  Section  [3.3.2], 

3.3.2   AMPX 

The  AMPX  Code  [G2]  is  a  modular  system  (NITAWL  and 
XSDRNPM  codes)  which  evaluates  123  group  cross-section  data 
with  resonance  parameters  to  generate  coupled  multi-group 
neutron  cross-sections.   The  reference  ball  design  used  for 
number  density  calculations  is  described  in  Appendix  [D.3]. 
The  reference  cell-averaged  ambient  temperature  was  900°C  [W3]. 

AI4PX  calculations  generated  all  cross-section  data  used 
in  the  present  work,  treating  the  46  nuclides  listed  in  Table 
[D.l].   For  the  depletion  chains  detailed  in  Table  [D.7],  appro- 
priate nuclide  decay  constants,  fission  product  yields,  energy 
releases  per  fission,  groupwise  x  spectra  and  the  resulting 
relative  groupwise  fluxes  are  given  in  Tables  [D.8],  [D.9], 
[D.IO],  [D.ll]  and  [D.12],  respectively,  as  functions  of  the 
carbon-to-heavy  metal  atom  ratio.   Energy  group  parameters  are 
given  in  Table  [D.13]. 

The  four-group  cross-section  data  generated  by  AMPX  and 
used  in  PDQ-7  calculations  is  listed  in  Appendix  C. 
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3.4   Computational  Methods  for  Evaluation  of  Homogeneous 
Zone  Parameters 

As  discussed  in  Chapter  2,  evaluation  of  breed/burn 

operation  in  a  homogeneous  one-zone  system  requires  computa- 

,  tion  of  system  criticality  and  fissile  enrichment  as  a  function 

! 

of  burnup.   The  performance  characteristics  of  an  arbitrary 
region  can  be  bracketed  through  the  evaluation  of  several 
reference  examples. 

Using  the  system  model  described  in  Section  [3.2],  the 
carbon- to-heavy  metal  atom  ratio  and  the  beginning-of-cycle 
fissile  enrichments  for  the  system  were  varied  in  a  series  of 
eigenvalue/depletion  calculations.   From  the  calculated 
information,  the  parameters  necessary  for  an  evaluation  of 
the  breed/burn  system  were  obtained. 

For  a  given  burn  zone  composition,  varying  the  initial 
fissile  enrichment  will  yield  a  series  of  similar  curves  when 
k   is  plotted  versus  the  fuel  burnup,  as  illustrated  in  Fig. 
[3.1a].   Achieving  a  value  of  k^  >  1.0  will  satisfy  the  first 
operating  criterion:  achieving  system  criticality.   Then,  a 
linear  burnup  model  of  the  form.: 


c         c 


can  be  applied  to  the  zone,  where  a  and  b  are  constants  and  n 
is  the  number  of  batch  fuel  loadings.   For  a  continuously- 
fueled  system,  whose  performance  can  be  approached  in  a  pebble 
bed  reactor,  n  ->  °°  so  that: 
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BURNUP ,  B 

FIGURE  3.1a;   REPRESENTATIVE  VARIATION  OF 
k   WITH  BURNUP  (linear) 


•» 


BURNUP ,  B 


FIGURE  3.1b:   REPRESENTATIVE  VARIATION  OF 
ENRICHMENT  WITH  BURNUP 
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^EOL^  =  2a(egQ^   -  b)  .  (3.2) 

c 


The  constants  a  and  b  can  be  calculated  if  batch  depletion 
calculations  are  performed  for  two  different  initial  enrich- 
ments; therefore,  to  achieve  any  design-limited  end-of-life 
burnup  for  continuous  burn  zone  refueling,  Eq .  (3.2)  can  be 
used  to  find  the  required  e^^^     .   Hence,  for  a  burn  zone  load- 

C 

ed  with  an  initial  fissile  enrichment  of  £,cv'  bred  in  a 
breed  zone,  the  maximum  obtainable  end-of -cycle  burnup  can  be 
calculated  using  Eq.  (3.1)  for  either  a  batch  load  or  for 
continuous  refueling. 

In  its  fullest  development,  the  breed/burn  cycle  uses 
the  heavy  metal  discharged  from  the  breed  zone  at  an  enrich- 
ment ^j^ay   3S  the  fuel  loaded  into  the  burn  zone  at  beginning- 
of-cycle,  again  at  an  enrichment  £,qv*   Therefore,  the  value 
of  £,CY  charged  to  the  burn  zone  must  be  high  enough  to  allov/ 
the  next  cycle  to  achieve  an  equivalent  breed  zone  discharge 
enrichment.   Several  successive  passes  through  the  breed  zone 
may  be  required  for  a  given  fuel  ball  to  achieve  a  desired 
enrichment,  but  as  a  very  rough  approximation,  assuming  one 
fissile  nucleus  is  bred  (and  preserved)  in  the  breed  zone  for 
every  fissile  nucleus  burned  in  the  burn  zone: 


bz  c 


^EOL,   '*>  i  ^ASY  '*'•  '3.4) 
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In  evaluating  this  criterion,  breed  zone  enrichments  at 
breed  zone  end-of-cycle  will  be  taken  to  be  equal  to  the 
asymptotic  values  given  in  Table  [2.5]  as  a  function  of  the 
ambient  spectrum. 

Burn  zone  enrichments  can  be  treated  more  explicitly. 
Because  the  PDQ-7  analysis  takes  into  account  the  effects  of 
U235  and  Pa233  production,  fission  product  buildup,  Th232 
depletion  and  parasitic  absorption,  all  of  which  were  neglect- 
ed in  the  simplistic  analysis  of  Chapter  2,  the  actual  system 
operating  characteristics  may  differ  slightly  from  those 
predicted  previously.   Concentrations  of  all  nuclides  appear- 
ing in  the  depletion  chains  are  calculated  by  PDQ-7  at  each 
time  increment.   Using  Eq.  (2.18),  the  enrichment  of  a  given 
burn  zone  as  a  function  of  burnup  can  be  calculated.   Operating 
at  a  constant  power  density  to  high  fuel  burnups  will  result 
in  fissile  enrichments  being  reduced  to  the  spectrum-averaged 
equilibrium  limits  at  end-of-cycle,  as  illustrated  in  Fig. 
[3.1b].   If  the  ambient  spectrum  in  the  breed  zone  is  much 
harder  than  that  in  the  burn  zone,  the  equilibrium  enrichment 
in  the  burn  zone  will  be  significantly  less  than  the  asymp- 
totic limit  achieved  in  the  breed  zone.   Burning  fuel  elements 
in  the  burn  zone  to  high  burnups  and  hence  to  low  fissile 
enrichments  at  end-of-life  (for  high  carbon- to-heavy  metal 
atom  ratio  compositions)  is  in  keeping  with  the  nonprolif era- 
tion  objectives  of  the  breed/burn  system.   Very  low  end-of-life 
fissile  enrichments  thus  allow  spent  fuel  elements  to  be 
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discharged  directly  to  waste  disposal,  in  accordance  with  the 
system  schematic  in  Fig.  [1.1], 

Zone  compositions  were  evaluated  using  PDQ-7  in  batch 
load  depletion  calculations,  so  the  maximum  burn  zone  burnups 
(achieved  using  continuous  refueling)  are  twice  that  obtained 
by  applying  Eq .  (3.1)  (setting  n  =  1)  to  batch  zone  calcula- 
tions . 

In  summary,  the  second  criterion  for  successful  operation, 
that  adequate  fissile  material  be  bred  in  the  breed  zone  to 
fuel  the  burn  zone,  is  accomplished  only  if  the  breed  zone 
end-of-life  enrichment  (here  estimated  as  ^j.^^)    is: 


bz  bz  c 


This  requirement,  combined  with  the  requirement  for  system 

criticality  are  necessary  (but  not  sufficient)  conditions 

for  successful  breed/burn  operation.   The  third  assessment 

criterion,  that  of  zone  leakage,  must  be  examined. 

Equation  (2.39)  defines  the  requirement  that  k^  >_   1.5. 

This  value  applies  only  when  the  burn  zone  breeding  ratio, 

b  ,  is  maximized,  which  occurs  under  end-of-life  asymptotic 

conditions,  or  when  b   =1.0.   Since  this  condition  (Eq. 

(2.39))  defines  a  minimum  value  for  k   ,  only  those  combina- 

tions  of  breed  and  burn  zone  compositions  for  which  the  first 

two  criteria  for  operation  have  been  met  and  for  which 

k      >  1.5  will  be  considered  as  feasible  possibilities. 
BOL 
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Should  all  three  criteria  be  met  by  any  given  candidate  system, 

then  a  more  detailed  analysis  examining  system  geometry  in  a 

^   two-zone  heterogeneous  configuration  will  be  required. 

For  this  analysis,  PDQ-7  depletion  calculations  were 

carried  out  for  beginning-of-lif e  enrichments  of  2.0%,  4.5% 

and  7.0%  U233,  for  zone  carbon-to-heavy  metal  atom  ratios  of 

110,  250,  325  and  450.   The  analyses  of  k   and  e,^„  as  a  func- 

^  °°       ASY 

tion  of  the  zone  burnup,  previously  described,  were  performed 
to  evaluate  the  operating  characteristics  of  the  breed/burn 
system.   Depletion  calculations  were  run  for  all  compositions 
through  fuel  burnups  of  120,000  MWD/MTHM.   A  summary  of  these 
results  is  given  in  Chapter  4 . 

3 . 5   Summary 

Procedures  used  to  analyze  the  characteristics  of  operat- 
ing breed/burn  fuel  cycles  in  pebble  bed  HTGRs  have  been  pre- 
sented.  A  series  of  PDQ-7  eigenvalue/depletion  runs  were 
carried  out  on  homogeneous  one-zone  systems  to  evaluate 
system  performance,  varying  the  parameters  of  the  carbon- to- 
heavy  metal  atom  ratio  (composition)  and  the  beginning-of-lif e 
fissile  enrichment. 

Successful  system  operation  will  be  achieved  only  if  three 
necessary  (but  not  sufficient)  conditions  are  met: 

(a)  achieving  system  criticality,  k^  >_  1.0, 

(b)  achieving  a  fissile  production  rate  such  that 


e.o^.      (^°)   >  B^^,      (%), 

"bz 
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and   (c)   accounting  for  zone  leakage  by  providing  adequate 
excess  reactivity,  k  >  1.5. 
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Satisfying  all  three  of  these  objectives  simultaneously  will 
allow  the  two-zone  layered  pebble  bed  system  to  operate  on  a 
stand-alone  breed/burn  fuel  cycle  (within  the  limits  of  the 
conservative  approximations  used  in  this  analysis) . 
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IV.   BREED/BURN  COMPUTATIONAL  RESULTS 

4 . 1  Introduction 

The  PDQ-7  runs  made  for  analysis  of  breed/burn  feasibility 
were  detailed  in  Chapter  3.   The  results,  to  be  described  in 
this  chapter,  will  characterize  the  operating  parameters, 
consistent  with  core  compositions  and  power  densities  as  a 
function  of  burnup  history  necessary  to  assess  the  capability 
of  the  pebble  bed  HTGR  system  to  achieve  breed/burn  operation. 
The  analysis  specifically  addresses  the  three  necessary  (but 
not  sufficient)  criteria  defined  for  successful  breed/burn 
operation: 

(a)  system  criticality,  k   >^  1.0, 

(b)  rate  of  fissile  production,  B      (%)  £  ^aSY    ^^^ 

c  bz 

and   (c)   accounting  for  zone  leakage,  k^  >^  1 . 5 . 

4.2  Variation  of  k^  with  Fuel  Burnup 

Three  depletion  cases  were  treated  for  each  homogeneous 
zone  composition  (carbon- to-heavy  metal  atom  ratios  of  110, 
250,  325  and  450),  corresponding  to  fissile  enrichments  at 
the  beginning-of-cycle  of  2.0%,  4.5%  and  7.0%.   For  each  case, 
the  depletion  time  increments  corresponded  to  fuel  burnups  of 
2000,  6000,  30,000,  60,000  90,000  and  120,000  MWD/MTHM.   Cur- 
rent materials  constraints  limit  cycle  burnups  to  an  upper 
value  of  roughly  120,000  MWD/MTHM  [Tl] . 
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For  the  twelve  cases  considered,  the  values  calculated 
for  k^  at  the  start  of  each  depletion  iteration  are  listed 
in  Appendix  B,  Table  [B.3].   As  illustrated  in  Figs.  [4.1- 
4.4],  it  can  be  seen  that  for  all  compositions,  as  expected, 
k^  decreases  with  fuel  burnup,  and  for  each  composition,  high- 
er initial  fissile  enrichment  results  in  higher  values  of 
k^  .   A  step  decrease  in  k^  is  noted  during  the  first  time 
increment  due  to  the  rapid  saturation  of  xenon-135.   All 
cases  exhibited  a  nearly  linear  relation  between  k   and 

CO 

burnup  beyond  60  00  MWD/MTHM. 

At  the  beginning  of  life  all  compositions  achieved  criti- 
cality,  with  the  exception  of  the  case  for  which  C/HM  =  110 
and  G^^-r  =  2.0%.   However,  all  cases  for  which  £„„^  =  2.0% 
were  subcritical  at  a  burnup  of  only  2000  MWD/MTHM  and  re- 
mained so  for  the  duration  of  the  cycle.   In  these  cases,  a 
rise  in  k   was  observed  at  mid-cycle,  apparently  due  to  an 
increase  in  fissile  content.   This  will  be  addressed  further 
in  the  next  section. 

The  linear  burnup  model  of  Eq.  (3.1)  can  be  applied  to 
each  burn  zone  case,  setting  n  =  1,  since  the  calculations 
were  done  for  batch  loads.   For  each  carbon- to-heavy  metal 
atom  ratio,  the  fissile  enrichments  necessary  at  the  beginning- 
of-cycle  necessary  to  achieve  any  given  end-of-cycle  burnup 
for  either  a  batch-loaded  or  continuously-fueled  system  can 
then  be  determined. 

A  linear  correlation  of  the  form: 
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FIGURE  4.1:   k^  VERSUS  FUEL  BURNUP  FOR  C/HM  =  110 
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FIGURE  4.2:   k   VERSUS  FUEL  BURNUP  FOR  C/HM  =  250 
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FIGURE  4.3:   k   VERSUS  FUEL  BURNUP  FOR  C/HM  =  325 
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FIGURE  4.4:   k   VERSUS  FUEL  BURN  FOR  C/HM  =  450 
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^oo   =  ^1  B  +  c^  (4.1) 

was  fit  (linear  least  squares)  to  the  linear  region  of  each 
curve  (i.e.  for  all  burnups  after  initial  Xel35  saturation)  and 
the  coefficients  c^  and  c.^   were  determined.   The  resulting  c 
and  c^  are  listed  in  Appendix  B,  Table  [B.5].   In  this  model, 
B  is  the  average  reactivity-limited  fuel  burnup. 

Equation  (4.1)  can  be  solved  for  the  burnup  at  which 
k   =  1.0,  B,  as  : 

CO  '       ' 

1.0  =  c-j^  B  +  c^   ,  (4.2) 

for  each  burn  zone  burnup  history  shown  in  Figs.  [4.1  -  4.4]. 
Since  initial  enrichments  of  2.0%  resulted  in  subcritical 
systems  (for  all  except  the  first  time  step)  for  all  carbon-to- 
heavy  metal  atom  ratios  examined,  these  cases  were  excluded 
from  further  analysis.   For  the  remaining  two  burnup  histories 
for  each  burn  zone  composition  (corresponding  to  initial  en- 
richments of  4.5%  and  7.0%)  the  linear  burnup  relation  was 
applied.   Simultaneous  solution  of  the  resulting  two  equations 
yields  the  coefficients  a  and  b  for  each  carbon- to-heavy  metal 
atom  ratio.   These  coefficients  are  listed  in  Appendix  B, 
Table  [B.6] . 

Thus,  for  any  initial  fissile  enrichment,  the  end-of- 
cycle  batch  fuel  burnups  can  be  calculated;  and  setting  n  ->■  «• 
in  Eq.  (3.1)  yields  the  burn  zone  burnup  limit  for  continuous 
refueling.   Table  [4.1]  lists  the  reactivity-limited  burnups 
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TABLE  4.1 
REACTIVITY-LIMITED  BURNUP*  FOR 
CONTINUOUSLY-FUELED  BURN  ZONES 

c       bz 

BURN  ZONE  C/HM  ATOM  RATIO 
110     250     325     450 


110 

BREED 

ZONE    250 

C/HM 

ATOM    325 

RATIO 

450 


67400  97590  94480  86950 

28210  67390  65092  58500 

19920  61000  58880  52480 

13140  55780  53790  47560 


*  Values  shown  are  MWD/MTHM 

t  Breed  Zone  e,qy  values  are  given  in  Table  [4.2] 


as  e 


EOC* 
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for  continuously-fueled  burn  zones  as  a  function  of  the  breed 

zone  asymptotic  limits  given  in  Table  [4.2].   These  asymptotic 

I 

'  enrichments,  being  the  average  end-of-cycle  enrichments,  e^^^^, 

EOC 

generated  in  the  PDQ-7  calculations  for  each  composition, 
approximate  the  asymptotic  values  predicted  for  each  spectrum, 

I 

as  can  be  seen  in  Table  [4.2] . 
I       The  burnups  obtained  vary  substantially  with  the  breed 
and  burn  carbon-to-heavy  metal  atom  ratio  combinations,  as 
shown  in  Fig.  [4.5].   The  highest  values  of  fuel  burnups  are 
obtained  for  a  breed  zone  for  C/HM  =  110  and  a  burn  zone 
composition  in  the  range  of  C/HM  =  250  to  450. 

23 

4 . 3   Variation  of  e        with  Fuel  Burnup 

Having  met  the  objective  of  achieving  system  criticality, 
the  next  objective  is  to  examine  the  fissile  production  poten- 
tial.  The  analysis  in  Chapter  2  indicated  that  a  maximum 
possible  asymptotic  fissile  enrichment  would  be  achieved, 
given  adequate  exposure  time,  dependent  on  the  ambient  spec- 
trum..  The  fissile  enrichments  (U233  only,  neglecting  U235) 
for  each  zone  composition  are  listed  as  a  function  of  fuel 
burnup  in  Table  [B.4].   As  illustrated  in  Figs.  [4.6  -  4.9], 
it  is  clear  that  an  asymptotic  limit  is  approached  at  very 
high  fuel  burnups,  independent  of  the  initial  zone  fissile 
enrichment.   A  comparison  is  made  in  Table  [4.2]  of  the  average 
U233  enrichment  at  end-of-cycle  inferred  from  these  results 
with  the  asymptotic  enrichment  limits  predicted  using  the 
simplistic  analysis  of  Chapter  2.   As  previously  discussed, 
predicted  values  of  £7,37  ^^^^^^  correspond  nearly,  but  not 
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TABLE  4 . 2 
COMPARISON  OF  ACTUAL  END-OF-CYCLE  U2  33 
ENRICHMENTS  AND  PREDICTED  ASYMPTOTIC  LIMITS 

C/HM  ATOM  RATIO 
110       250       325       450 


z^*         0-03526   0.02686   0.02497   0.02302 


e        0.03697   0.02137   0.01807   0.01537 

^ASY 

iTT"     1.048     0.7956    0.7237    0.6677 


*  Predicted  using 


2  3 

a 
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PEAK  BURNUP  REGION 


FIGURE  4.5:   REACTIVITY  LIMITED  BURNUP  FOR 
ASYMPTOTIC  ENRICHMENT  LOADINGS 
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FIGURE  4.6:  U233  ENRICHMENT  VERSUS  FUEL  BURNUP 
FOR  C/HM  =  110 
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FIGURE  4.7:   U233  ENRICHMENT  VERSUS  FUEL  BURNUP 
FOR  C/HM  =2  50 
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FIGURE  4.8:   U233  ENRICHMENT  VERSUS  FUEL  BURNUP 
FOR  C/HM  =  3  25 
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FIGURE  4.9:   U233  ENRICHMENT  VERSUS  FUEL  BURNUP 
for  C/HM  =  4  50 


-91- 


exactly,  to  £gQj^/  due  to  the  more  rigorous  evaluation  of 
other  concurrent  effects  by  PDQ-7. 

A  breed  zone,  initially  having  no  fissile  content,  would 
eventually  build  up  to  a  fissile  concentration  near  the 
asymptotic  limit,  as  indicated  by  the  low-enriched  (2.0%  U233) 
depletion  cases  (which  also  exhibited  an  increase  in  k   with 
burnup  due  to  the  high  fissile  production  rate) .   It  was  noted 
in  the  previous  section  that  an  optimum  combination  of  zone 
compositions  to  maximize  obtainable  burn  zone  burnup  would 
use  a  breed  zone  where  C/HM  =  110  and  a  burn  zone  where  C/HM 
=  250  to  450.   Continuous  refueling,  which  can  be  approached 
in  pebble  bed  systems,  would  allow  burn  zone  burnups  in  these 
cases  approaching  100  GWD/MTHM.   Recognizing  that  for  contin- 
uous fueling  very  high  burnups  could  be  attained  (Table  [4.2]), 

the  criteria  for  breed  zone  burnup,  z^^^       <   B        ,  can  be 

c       bz 
achieved.   Additionally,  it  is  clear  from  Appendix  B,  Table 

[B.3],  that  increasing  the  carbon- to-heavy  metal  atom  ratio 

for  a  given  initial  fissile  enrichment  will  increase  k^ 

BOL 
This  effect  is  illustrated  in  Fig.  [4.10].   Because  system 

criticality  must  be  maintained,  higher  burn  zone  carbon-to- 
heavy  metal  atom  ratios  (^400  or  greater)  would  therefore 
be  preferred,  even  though  the  obtainable  burn  zone  burnups 
are  not  quite  the  maximum  possible. 

Breed/burn  system  design  in  the  pebble  bed  HTGR  can, 
with  careful  selection  of  zone  compositions,  meet  the  rather 
conservative  criterion  established  here  for  potential  fissile 
production  and  hence  breed/burn  cycle  operation  may  be 
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FIGURE  4.10:   k      VERSUS  C/HM  ATOM  RATIO  AS  A 
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feasible.   In  this  analysis,  the  effects  of  poisons,  control 
materials  and  reactor  vessel  materials  were  neglected  by 
analyzing  an  infinite  homogeneous  system.   Therefore,  these 
results  are,  in  fact,  extremely  conservative  estimates,  indi- 
cating that  true  system  performance  might  be  substantially 
poorer.   The  operation  of  a  two-zone  system,  however,  as 
described  in  Chapter  3,  requires  an  evaluation  which  accounts 
for  zone- to- zone  leakage. 

4 . 4   Criticality  with  Zone  Leakage 

Equation  (2.39)  required,  as  a  minimum,  that  at  the 
beginning-of-lif e: 


k^   >  1.5  .  (2.39) 

c 


This  was  necessary  to  account  for  the  significant  neutron 
leakage  out  of  the  burn  zone  (into  the  breed  zone) .   For  the 
reference  homogeneous  systems  evaluated,  beginning-of-lif e 
values  of  k   exceeded  1.5  only  for  initial  enrichments  of  7.0%. 
Since  the  asymptotic  limit  for  fissile  enrichment  in  a  breed 
zone  for  which  C/HM  =  110  is   3.5%,  as  shown  in  Fig.  [4.10], 
it  can  be  seen  that  idealized  core  configurations  yield  no 
values  of  k      >  1.5.   This  indicates  that  it  is  not  possible 

oo        ^ 

BOL 
to  achieve  a  k   high  enough  to  accomodate  zone  leakage  when 

using  breed  zone  discharge  fuel  in  the  burn  zone  (depicted  as 

the  "feasible"  range  in  Fig.  [4.10]).   Thus,  no  combination 

of  the  zone  compositions  evaluated  in  the  present  work  are 
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capable  of  satisfying  all  three  criteria  for  successful  breed/ 
burn  operation.   For  these  compositions,  then,  it  must  be 
concluded  that  it  is  not  feasible  to  operate  a  pebble  bed 
HTGR  system  on  the  breed/burn  fuel  cycle.   Moreover,  since  k 

-■  'CO 

appears  to  attain  its  maximum  value  (for  all  enrichments) 
at  a  carbon-to-heavy  metal  C/HM  value  of  '^^  4  50,  it  does  not 
appear  to  be  useful  to  consider  the  use  of  even  softer  spectrum 
burn  zones.   Referring  back  to  Table  [2.5],  extrapolation  to 
C/HM  <  110  indicates  that  a  harder  spectrum  breed  zone,  if 
physically  practicable,  might  have  an  asymptotic  enrichment 
greater  than  the  5%  or  so,  which  would  permit  k   >  1.5  when 
cycled  into  the  burn  zone. 

4 . 5   Summary 

Analysis  of  the  results  of  PDO-7  depletion  calculations 
for  homogeneous  zones  has  indicated  that  with  high  fuel  burnups, 
breed  zone  fissile  enrichments  approach  steady  state  values, 
as  predicted  in  the  analysis  of  Chapter  2.   These  enrichments, 
when  appropriate  zone  compositions  are  employed,  will  allow 
burn  zones  of  infinite  extent  to  achieve  adequately  high  burn- 
ups while  maintaining  k   >^  1.0.   Considering  only  these  cri- 
teria, breed/burn  system  operation  appears  feasible.   However, 
requiring  the  system  to  overcome  the  effects  of  burn  zone  leak- 
age showed  that  none  of  the  breed  and  burn  zone  combinations 
examined  in  the  present  work  could  achieve  adequate  values  of 
k   .   Hence,  the  potential  for  successful  breed/burn  operation 
in  pebble  bed  reactors  in  not  confirmed. 
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j       It  may  be  possible,  however,  by  further  reductions  in  the 
I  carbon- to-heavy  metal  ratio,  in  the  breed  zone,  to  increase 
the  asymptotic  enrichment  limit  of  the  breed  zone  (to  >  5%, 
I  and  preferably  ^   7%)  and  thereby  make  it  possible  to  achieve 

a  system  k^  ^  1.5.   Additionally,  as  can  be  seen  in  Fig.  [4.4], 
■   the  use  of  slightly  higher  carbon- to-heavy  metal  ratios  in 
j   the  burn  zone  (say  %  500)  can  produce  slightly  higher  values 
of  k^  ^    and  hence  it  may  be  possible  to  satisfy  the  leakage 
criterion.   Reasonable  compositions  must  be  maintained,  hov;- 
ever,  because  as  the  carbon-to-heavy  metal  atom  ratio  increases, 
the  required  core  volume  will  increase  for  a  given  thermal 
output  requirement. 

Another  option  would  be  to  utilize  the  pebble  bed  HTGR 
in  a  near-breed/burn  mode,  by  providing  the  additional  fissile 
makeup  necessary  to  achieve  the  required  values  of  k   .   Reac- 
tivity-limited burnups  would  thereby  be  increased,  thus  allow- 
ing the  system  to  better  utilize  its  potential  for  U233  pro- 
duction. 

The  goals  of  non-proliferation,  while  best-achieved  by 
a  successful  "stand-alone"  breed/burn  system,  might  still 
have  acceptable  operating  parameters  for  a  near-breed/burn 
system.   Fissile  U235  or  U233  loaded  into  the  system  in  a 
denatured  form,  would  not  violate  non-proliferation  standards 
and  could  provide  the  necessary  reactivity  to  effect  success- 
ful breed/burn  operation  in  the  pebble  bed  HTGR  system.   Should 
such  a  modified  breed/burn  fuel  cycle  be  considered,  however, 
the  effects  of  plutonium  production  would  be  to  be  evaluated. 
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V.   SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

5 . 1   Introduction 

The  realization  that  a  growing  number  of  nonnuclear 
weapons  states,  desiring  to  provide  for  their  own  growing 
energy  needs,  have  opted  to  develop  civilian  nuclear  power 
programs,  and  the  implied  relationship  between  nuclear 
electric  plants  and  weapons,  has  resulted  in  a  broader  inter- 
pretation of  the  implications  of  proliferation.   Both  national 
and  international  reassessments,  e.g.  the  NASAP  and  INFCE 
efforts,  are  underway  to  prepare  new  guidelines  for  the  devel- 
opment and  application  of  the  nuclear  fuel  cycle  for  commercial 
electrical  power  production. 

The  transition  to  a  breeder  economy  is  not  in  prospect  in 
the  United  States,  under  current  governmental  policy,  and 
alternative  options  must  utilize  nuclear  fuel  in  a  manner  so 
as  to  minimize  the  production  of  and/or  commercial  traffic  in 
materials  suitable  for  weapons  applications.   The  present 
thesis  is  concerned  with  one  such  class  of  proliferation- 
resistant  fuel  cycles,  designated  the  "breed/burn"  concept. 

Fischer  et  al  [Fl]  proposed  the  implementation  of  a  breed/ 
burn  system  in  a  fast  reactor  core,  using  depleted  uranium 
monocarbide  feed  fuel.   After  breeding  sufficient  fissile 
material  in  a  hard  spectrum  "breed  zone,"  fuel  assemblies  are 
moved  to  a  more  highly  moderated  region  of  the  core  where 
energy  is  generated  by  the  burning  of  the  now  slightly 
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enriched  fuel. 

Application  of  this  concept  to  thermal/epithermal  reactor 
systems  requires  the  use  of  a  compositionally  versatile  system, 
capable  of  achieving  the  high  fuel  burnups  necessary  for  effec- 
tive fissile  conversion. 

The  pebble  bed  HTGR  offers  inherent  versatility  in  fueling 
strategies,  core  burnup,  fuel  design  and  neutron  spectrum 
tailoring,  as  well  as  the  ability  to  approach  continuous  fuel- 
ing operation  by  using  flow- through  schemes.   Hence,  the  pre- 
sent work  focuses  primarily  on  the  utilization  of  the  U233/Th 
fuel  cycle  in  a  breed/burn  configuration  of  the  pebble  bed 
HTGR. 

5 . 2   Summary  of  Analysis 

Fischer  et  al  [Fl]  have  proposed  that  in  implementing  a 
fast  reactor  breed/burn  fuel  cycle,  the  primary  objective  is 
to  maximize  the  breeding  ratio  through  spectrum  hardening. 
More  highly-moderated  regions  are  used  for  fuel  burning,  but 
a  hardened  ambient  spectrum  in  the  breed  zone  is  required, 
presumably  to  increase  the  overall  system  breeding  ratio. 
Primarily  through  decreases  in  parasitic  absorption  (in 
fissile  as  well  as  structural  material  and  fissile  products) 
and  increases  in  fertile  fission  rates,  the  breeding  ratio 
does  increase  with  hardening,  but  it  was  found  in  the  present 
work  that  this  is  accompanied  by  a  decrease  in  the  asymptotic 
fissile  enrichment  obtained  for  ultra-high  burnups,  when  the 
rate  of  fertile  capture  equals  the  fissile  absorption  rate. 
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This  spectrum-dependent  apparent  limit  is  expressed  by  the 
relation: 

,   fertile,   fissile> 

C  Pi 

^ASY  "  I   +    f^g    fertile  ,^  fissile  ^^'^^ 

For  fast  reactor  applications,  the  preference  is  for 
softer  breed  zone  spectra:  softer  than  a  typical  LMFBR  but 
still  harder  than  in  common  thermal  reactors. 

Developing  the  breed/burn  fuel  cycle  for  thermal  reactor 
configurations  necessitates  the  use  of  a  system  capable  of 
achieving  regionwise  spectral  shifts  (in  order  to  segregate 
breed  and  burn  zones)  while  being  able  to  achieve  high  conver- 
sion ratios.   The  HTGR  is  the  most  versatile  among  thermal/ 
epithermal  reactors  currently  of  commercial  interest.   This 
inherently  superior  system,  when  employed  with  the  U233/Th 
fuel  cycle  (chosen  because  of  its  outstanding  neutronic 
properties  in  the  thermal  and  epithermal  energy  ranges) 
affords  the  required  flexibility  for  implementation  of  a  breed/ 
burn  cycle.   Very  high  fuel  burnups,  relatively  low  enrichment 
requirements  and  the  stated  requirements  for  zoning  flexibility 
made  the  pebble  bed  HTGR  a  logical  system  choice.   The  present 
work  focused  primarily  on  the  use  of  this  system  for  obtaining 
a  feasible  thermal/epithermal  breed/burn  design. 

Contrary  to  the  situation  with  fast  reactor  breed/burn 
systems,  increasing  the  breed  region  asymptotic  fissile  en- 
richment in  the  U233/Th  pebble  bed  HTGR  is  accomplished  by 
spectrum  hardening.   In  addition,  the  k   values  for  a  given 
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enrichment  increase  as  the  moderator-to-fuel  ratio  (carbon- 
to-heavy  metal)  is  increased.   Therefore,  a  system  utilizing 
two  zones;  a  very  hard  breed  zone  and  a  very  thermal  burn 
region  appears  most  promising. 

Successful  operation  of  any  breed/burn  system  requires 
that:   (1)  system  criticality  be  maintained,  or  k   >   1.0, 

CO      — 

and  (2)  breed  zone  fissile  production  be  adequate  to  fuel 

the  burn  zone,  or,  using  a  very  rough  approximation  for  the 

breed  zone  discharge  enrichment,  e      (%)  ~  £,,„„   (%)  = 

bz  bx 

B      (%) .   Burn  zone  operation  must  then  require  an  initial 


c 


enrichment,  ^dqt  /  J^o  greater  than  that  discharged  from  the 

c 
breed  zone,  or  eyvgy    (%)  t   ^fot   ^^^  *  ^   final  operating 

bz  c 

criterion  requires  that  system  criticality  consider  the  effects 

of  zone-to-zone  neutron  leakage.  In  applying  a  leakage  correc- 
tion to  the  macroscopic  absorption  cross-section,  the  criterion 
for  criticality  yields  the  necessity  that: 


-  1  >  1  +  ^   ,  (5.2) 

a 


where  e   is  the  "leakage  cross-section."   By  further  manipula- 
tion  of  Eq.  (5.2),  using  conservative  (optimistic)  approxima- 
tions, one  finds  that: 


i 


k    -  1  >  .  /,  (5.3) 

oo       —  1  +  b 
c 


and  since  the  burn  zone  breeding  ratio,  b^,  reaches  a  maximum 
of  unity  at  the  asymptotic  limit,  for  successful  operation, 
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I  accounting  for  zone-to-zone  leakage^ k^  must  be  greater  than 
1.5. 

Three  operating  criteria  have  thus  been  established  for 
successful  breed/burn  operation.   These  necessary  (but  not 
sufficient)  conditions  define  the  feasible  range  of  operation 
for  the  system. 

To  fully  evaluate  the  practicality  of  pebble  bed  systems, 
state-of-the-art  computer  calculations  were  performed  using 
the  PDQ-7  code  [CI] .   Four-group  cross-section  data  was  gen- 
erated for  the  pebble  bed  system  using  the  AMPX  modular  code 
system  (NITAWL  and  XSDRNPM  codes) .   Homogeneous  one- zone 
regions  were  evaluated  in  a  series  of  eigenvalue/depletion 
calculations,  varying  the  zone  composition  (carbon-to-heavy 
metal  atom  ratios)  and  initial  enrichments.   Depletion  runs 
were  performed  for  fuel  burnups  to  120,000  MWD/MTHM. 

For  the  various  compositions  evaluated  in  this  thesis, 
asymptotic  zone  enrichments  are  shown  in  the  first  column  of 
Table  [5.1].   Increasing  asymptotic  enrichments  are  observed 
for  decreasing  carbon- to-heavy  metal  atom  ratios,  justifying 
the  preference  for  a  hardened  breed  zone  spectrum  in  this 
reactor  type.   Using  the  highest  available  asymptotic  enrich- 


ment, e^sY 


,  as  the  initial  fissile  enrichment 


C/HM  =  110 
in  a  continuously-refueled  burn  zone,  the  burn  zone  end-of- 

life  burnups  are  also  shown  in  the  first  row  of  Table  [5.1]^ 

as  a  function  of  the  zone  carbon- to-heavy  metal  atom  ratio. 

It  is  clear  that  this  reactor  system  is  able  to  achieve  high 
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TABLE  5.1 

PERFORMANCE  PARAMETERS  AS  A  FUNCTION  OF 
BREED  AND  BURN  ZONE  COMPOSITION 


BURN  ZONE  C/HM  ATOM  RATIO 
110        250        325        450 


B 


EOL  , 
c' 


'ASY,   , 


6.740 


9.759 


k+ 


BOL 


9.448 


8.695 


BREED 

ZONE 

C/HM 

ATOM 

RATIO 


110 


250 


325 


450 


0.0353 
(0.0370) 

0.0269 
(0.0214) 

0.0250 
(0.0181) 

0.0230 
(0.0154) 


1.069 
(1.110) 

0.971 
(0.952) 

0.948 
(0.910) 

0.924 
(0.865) 


1.246 
(1.256) 

1.134 
(1.095) 

1.107 
(1.052) 

1.080 
(1.005) 


1.290 
(1.286) 

1.171 
(1.125) 

1.143 
(1.082) 

1.113 
(1.035) 


1.323 
(1.313) 

1.196 
(1.152) 

1.166 
(1.109) 

1.135 
(1.062) 


^^^  MWD/MTHM 


In 


,  using  e 


10 


BOL 


ASY 


bz  C/HM  =  110 


continuous  refueling. 


t 


for 


For  comparison,  numbers  in  parenthesis  were  predicted  by  the 
one-group  calculation  of  Chapter  2. 
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infinite  medium  fuel  burnups  while  requiring  quite  low  begin- 
ning-of-life  enrichments. 

Values  for  k^     are  given  in  Table  [5.1]  as  a  function 

BOL 
c 

of  the  combination  of  breed  and  burn  zone  compositions.   Al- 
though most  combinations  result  in  criticality  for  the  infinite 
burn  zone  (k^>1.0),   none  of  the  combinations  of  zone  composition 
reach  k^   values  of  1.5.   Hence,  based  on  this  simplistic 
analysis,  none  of  the  compositions  evaluated  in  this  thesis 
are  capable  of  overcoming  the  zone-to-zone  leakage  effects. 

5 . 3   Conclusions 

Successful  stand-alone  breed/burn  operation  could  not 
be  accomplished  in  U233/Th  pebble  bed  HTGR  systems  for  the 
ranges  of  compositions  considered  in  this  analysis: 
110  <_   C/HM  £  450.   While  the  pebble  bed  system  provides  the 
flexibility  needed  to  implement  a  successful  fuel  cycle,  the 
effects  of  zone-to-zone  leakage  necessitates  values  of  k 
that  are  too  high  to  be  achieved  with  breed  zone  asymptotic 
enrichments  (<3.5%  in  the  C/HM  range  studied).   Therefore,  if 
the  HTGR  is  still  to  be  considered  as  a  viable  breed/burn 
system,  additional  breed  zone  hardening  will  be  necessary  to 
increase  the  attainable  fisssile  enrichment:  carbon- to-heavy 
metal  atom  ratios  much  less  than  110  must  be  contemplated,  so 
as  to  achieve  asymptotic  enrichments  on  the  order  or  7%  or  so. 
Increasing  the  carbon- to-heavy  metal  atom  ratio  in  the  burn 
zone  could  possibly  be  used  to  increase  k_^  ,  thereby  helping 
to  overcome  the  constraints  imposed  by  the  zone-to- zone 
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leakage.   Such  an  action  must  consider,  however,  the  additional 
reactor  volume  required  for  a  given  thermal  output,  and  the 
fact  that  k^  approaches  a  broad  maximum  just  beyond  C/HM-450. 

For  fast  reactor  applications,  the  breed/burn  concept 
offers  good  potential  for  successful  operation.   However, 
"soft"  fast  reactor  spectra  would  be  most  likely  to  meet  all 
operating  criteria,  contrary  to  the  proposal  by  Fischer  et  al 
[Fl] .   The  much  higher  asymptotic  enrichments  available  in 
these  spectra  could  allow  a  breed/burn  system  to  produce  large 
fissile  masses  in  the  breed  zone  in  a  time  span  consistent 
with  lower  fuel  burnups  in  the  burn  zone. 

A  successful  breed/burn  system,  can  satisfy  the  nonproli- 
feration  objectives  of  minimizing  the  traffic  in  sensitive 
materials  suitable  for  weapons  application  at  all  points  in 
the  fuel  cycle.   Additionally,  the  use  of  fertile-only  makeup 
makes  fuel  enrichment  unnecessary,  while  the  high  burn  zone 
burnups  attainable  help  reduce  discharge  fissile  enrichments 
to  levels  which  are  insufficient  to  justify  reprocessing. 

5 . 4   Recommendations  for  Future  Work 

The  results  of  this  analysis  showed  asymptotic  enrich- 
ments in  fast  reactor  spectra  were  at  their  highest  levels 
in  the  "soft"  end  of  the  available  spectra.   On  the  other 
hand,  the  highest  asymptotic  enrichments  in  the  epi thermal 
range  were  obtained  in  the  hardest  available  epithermal  spec- 
tra.  An  apparent  maximum  exists,  therefore,  for  a  reactor 
operating  in  the  range  between  a  "soft"  LMFBR  and  a  "hard" 
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HTGR.   Future  analyses  should  examine  the  possibility  of  de- 
signing a  system  to  operate  in  the  described  range,  thereby 
maximizing  the  potential  for  successful  breed/burn  operation. 

In  addition,  Karam  [Kl]  has  proposed  the  use  of  suspended- 
bed  reactors  (SBR)  as  alternatives  to  LMFBRs .   These  reactors 
would  utilize  pyrolytic  carbon-coated  dicarbide  or  metallic 
fuel  kernels  (i.e.  HTGR  fuel  particles)  suspended  by  upward 
flowing  helium  gas  (coolant)  in  a  cluster  of  pressurized  tubes. 
High  breeding  ratios  (1.2  to  1.38)  are  projected  for  this 
system,  since  the  utilization  of  HTGR-type  fuel  kernels  v;ithout 
the  surrounding  graphite  matrix  results  in  an  extremely  hard 
epi thermal  spectrum.   Thus,  it  may  be  possible  to  adapt  some 
of  the  features  of  the  SBR  to  achieve  a  spectrum  in  the  range 
required  for  maximizing  the  asymptotic  enrichment  limit.   Com- 
bining the  high  breeding  ratios  and  hard  spectrum,  and  recog- 
nizing the  ability  to  maintain  composition  zoning  through  the 
use  of  tube-contained  fuel,  this  type  system  should  be  investi- 
gated carefully  for  its  applicability  to  the  breed/burn  system 
concept. 

Less  exotic  HTGRs  still  offer  exceptional  potential  for 
the  implementation  of  alternate  fueling  schemes.   Even  apart 
from  the  SBR,  conventional  or  pebble  bed  type  HTGR  composi- 
tions for  which  C/HM  <  110  should  be  investigated,  as  it  may 
be  possible  to  increase  the  available  asymptotic  enrichments 
to  levels  high  enough  to  produce  the  necessary  burn  zone  excess 
reactivity  to  overcome  zone-to-zone  leakage  effects.   Only 
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after  thorough  analysis  of  all  possible  compositions  should 
the  HTGR  be  judged  as  either  feasible  or  unfeasible  for  imple- 
mentation of  a  breed/burn  fuel  cycle. 

Finally,  once  a  combination  of  breed/burn  zones  pass  the 
preliminary  requirements  proposed  in  this  thesis,  analyses 
should  proceed  to  the  logical  step  of  multi-zone  PDQ-7  calcu- 
lations in  which  the  zone-to-zone  coupling  is  explicitly 
accounted  for.   This  is  essential  if  the  extent  to  which  zones 
can  be  spectrally  isolated  is  to  be  established  --  an  essential 
prerequisite  for  successful  breed/burn  operation. 
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APPENDIX  A 
ONE-GROUP  CONSTANTS 

In  this  appendix  are  tabulated  one-group  cross-sections 
(Tables  [A. 2]  through  [A. 5])  collapsed  from  the  four-group 
cross-sections  given  in  Appendix  C.   The  four-group  cross- 
sections,  generated  using  the  AMPX  Code  System  [G2] ,  were 
flux  weighted  using  the  relative  groupwise  flux  values 
listed  in  table  [A.l].   Nuclides  are  identified  in  the  cross- 
section  tables  by  an  identification  number,  "I#",  listed 
alphabetically  by  isotope  in  Table  [D.l]. 
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TABLE  A.l 
RELATIVE  GROUPWISE  NEUTRON  FLUX* 

C/HM  ATOM  RATIO 
GROUP    Au      110     250     325     450 

1  4.400      1.000      1.000      1.000      1.000 

2  5.750      1.534      1.533      1.533      1.534 

3  5.750      1.025      1.194      1.236      1.282 

4  5.916      0.920      2.305      2.939      3.874 


*  Normalized  to  ({)-,  for  each  carbon-to-heavy 
metal  atom  ratio. 
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APPENDIX  B 
TABULATION  OF  DATA  USED  IN  CHAPTERS  2  AND  4 

In  this  appendix,  data  is  tabulated  for  the  Sheaffer 
et  al  [Si]  cross-section  correlation  (Tables  [B.l]  and  [B.2]) 
used  in  Chapter  2,    as  well  as  the  data  resulting  from  the 
computer  runs  specified  in  Chapter  3  and  discussed  in  Chapter 
4,  (Tables  [B.3]  through  [B.6]). 
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TABLE  B. 2 
NEUTRON  CROSS-SECTIONS  FOR  FAST  SPECTRA* 


s 

a 

<' 

°l' 

°r 

<' 

°£' 

0.3 

0.6028 

3.739 

3.214 

0.4556 

2.623 

2.050 

0.4 

0.4550 

3.306 

2.907 

0.3806 

2.401 

1.964 

0.5 

0.3658 

3.005 

2.689 

0.3311 

2.242 

1.900 

0.6 

0.3061 

2.780 

2.523 

0.2954 

2.121 

1.849 

0.7 

0.2633 

2.603 

2.391 

0.2683 

2.023 

1.807 

0.8 

0.2311 

2.459 

2.282 

0.2468 

1.941 

1.771 

0.9 

0.2060 

2.338 

2.190 

0.2217 

1.872 

1.740 

1.0 

0.1858 

2.235 

2.111 

0.2147 

1.813 

1.713 

*From  Sheaf fer  et  al  [SI]  correlation,  o   =  o ^S^ . 
Cross-sections  in  barns. 
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TABLE  B.3 

k^  AS  A  FUNCTION  OF  FUEL  BURNUP  (from  PDQ-7) 

C/HM  ATOM  RATIO 

110  250      325  450 

0.8863  1.036  1.067  1.086 

0.8226  0.9450  0.9705  0.9853 

0.8356  0.9038  0.9120  0.9072 

0.9114  0.9091  0.8837  0.8335 

0.9169  0.8949  0.8713  0.8343 

0.9046  0.8724  0.8477  0.8102 

1.164  1.353  1.402  1.442 

1.103  1.269  1.313  1.350 

1.081  1.231  1.270  1.302 

1.010  1.086  1.095  1.090 

0.9529  0.9515  0.9288  0.8851 

0.9180  0.8827  0.8523  0.8062 

1.307  1.495  1.548  1.597 

1.253  1.416  1.463  1.506 

1.230  1.386  1.431  1.473 

1.133  1.257  1.290  1.315 

1.036  1.098  1.103  1.094 

0.9651  0.9601  0.9352  0.8893 


BOL 

BURNUP 

(%) 

(MWD/MTHM) 

2.0 

0 

2000 

6000 

30,000 

60,000 

90,000 

4.5 

0 

2000 

6000 

30,000 

60,000 

90,000 

7.0 

0 

2000 

6000 

30,000 

60,000 

90,000 

I 


-116- 


TABLE  B.  4 
U2  3  3  ENRICHMENT  AS  A  FUNCTION  OF  FUEL  BURNUP  (from  PDQ-7) 


BOL 

BURNUP 

C/HM  ATOM  RATIO 

(%) 

(MWD/MTHM) 

110 

250 

325 

450 

2.0 

0 

2.0 

2.0 

2.0 

2.0 

2000 

1.937 

1.845 

1.828 

1.814 

6000 

2.114 

1.758 

1.681 

1.603 

30,000 

3.047 

2.019 

1.771 

1.521 

60,000 

3.422 

2.116 

1.846 

1.615 

90,000 

3.524 

2.095 

1.813 

1.574 

120,000 

3.513 

2.038 

1.753 

1.512 

4.5 

0 

4.5 

4.5 

4.5 

4.5 

2000 

4.367 

4.310 

4.300 

4.292 

6000 

4.291 

4.061 

4.013 

3.964 

30,000 

4.044 

3.168 

2.944 

2.711 

60,000 

3.844 

2.487 

2.158 

1.845 

90,000 

3.720 

2.192 

1.865 

1.581 

120,000 

3.616 

2.072 

1.772 

1.524 

7.0 

0 

7.0 

7.0 

7.0 

7.0 

2000 

6.840 

6.799 

6.792 

6.785 

6000 

6.667 

6.494 

6.458 

6.419 

30,000 

5.810 

5.115 

4.929 

4.729 

60,000 

4.956 

3.719 

3.388 

3.059 

90,000 

4.358 

2.772 

2.370 

1.990 

120,000 

3.962 

2.272 

1.897 

1.575 
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TABLE  B.5 
CURVE-FIT  COEFFICIENTS  FOR  LINEAR 
k   VERSUS  BURNUP  CORRELATION 


e„^^   COEFFICIENT  C/HM  ATOM  RATIO 

BOL 

(%)  110        250        325        450 

4.5        c  1.092      1.251  1.292  1.326 

c^  -2.092E-6  -4.445E-6  -5.336E-6  -6.359E-6 

7.0        Cq  i;247      1.418  1.470  1.520 

c,  -3.286E-6  -5.177E-6  -5.995E-6  -7.009E-6 


k   =  c,B  +  Cp, 

oo       1         U 
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TABLE  B . 6 
CURVE-FIT  COEFFICIENTS  FOR 
LINEAR  BURNUP  MODEL 

COEFFICIENT  C/HM  ATOM  RATIO 

110  250        325        450 

a        1.256E+6  9.680E+5   9.420E+5   9.118E+5 

b        1.014E-2  -1.344E-2  -1.318E-2  -1.071E-2 


B^^^  =  a{e^^^    -  b) 


BOL     ^  BOL     ' n  +  1 
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APPENDIX  C 
FOUR-GROUP  CONSTANTS 

In  this  appendix  are  tabulated  the  four-group  cross- 
section  data  generated  by  the  AMPX  code  system  for  carbon- 
to-heavy  metal  atom  ratio  of  110,  250,  325,  and  450.   Refer 
to  Appendix  D  and  Section  [3.3.2]  for  code  descriptions  and 
input  parameters.   Table  [D.13]  shows  the  neutron  energy 
group  structure  and  Table  [D.l]  gives  the  identity  of  each 
nuclide  corresponding  to  the  I#  in  the  tables  of  this 
appendix. 
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APPENDIX  D 
COMPUTER  CODE  INPUT  PARAMETERS 

D. 1   Introduction 

In  this  appendix  are  detailed  the  parameters  used  in  the 
PDQ-7  calculations  described  in  Chapter  3.   Specific  spectrum- 
dependent  parameters  and  nuclide  depletion  chain  paramters 
were  generated  using  the  AMPX  code. 

D. 2   Nuclide  Identification 

Table  [D.l]  lists  the  46  nuclides  treated  by  the  AMPX 
code  system.   Reference  numbers  are  used  for  identification  of 
nuclides  in  PDQ-7. 

D. 3   Number  Densities 

This  section  describes  the  reference  ball  design  used 
both  in  cross-section  generation  in  the  AMPX  code  and  in 
modeling  the  pebble  bed  breed/burn  system  for  PDQ-7. 
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TABLE  D.l 
NUCLIDE  IDENTIFICATION 


NUCLIDE 

REFERENCE 
NUMBER^ 

I#^ 

NUCLIDE 

REFERENCE 
NUMBER^ 

i#5 

. Am243 

24395 

20 

Pml47 

14761 

12 

Carbon 

1206 

2 

Pinl4  8 

14861 

14 

Csl33 

13355 

9 

Pinl48m 

14800 

13 

Csl34 

13455 

10 

Pml49 

14461 

28 

Csl35 

13555 

11 

Pa233 
Recvel 

23391 

42 

Eul53 

15363 

17 

1200 

6 

Eul54 

15463 

40 

Rhl03 

10345 

31 

-,Eul55 

15563 

41 

Rhl0  5 

10545 

32 

^FP  1 

10001 

22 

Sml49 

14962 

15 

FP  2 

10002 

21 

Sml50 

15062 

38 

Helium 

402 

1 

Sinl51 

15162 

16 

1135 

13553 

4 

Sml52 

15262 

39 

Kr82 

8236 

29 

Silicon 

1428 

7 

Kr83 

8336 

30 

Agl09 

10947 

33 

Ndl43 

14360 

36 

Th232 

23290 

23 

Ndl45 

14560 

37 

U233 

23392 

24 

Ndl47 

14760 

3 

U234 

23492 

25 

Np239 

23993 

27 

U235 

23592 

26 

Oxygen 

1608 

5 

U236 

23692 

18 

Pu239 

23994 

44 

U238 

23892 

43 

Pu240 

24094 

45 

Xel31 

13154 

34 

Pu241 

24194 

46 

Xel33 

13354 

35 

Pu242 

24294 

19 

Xel35 

13554 

8 

Slowly-saturating  fission  products  (lumped) 

Non-saturating  fission  products  (lumped) 

Reciprocal-velocity  correction  to  account  for  graphite 
impurities  [W3] 

Combined  mass  number  and  atomic  number  for  identification 
in  PDQ-7 

'identification  number  for  AMPX  code  and  cross-section 
tables  of  Appendices  [A]  and  [C] 
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D.3.1   Nomenclature 


C/HM  =  Ratio  of  carbon  number  density  to  heavy  metal 
number  density 

e     =  Fissile  enrichment 


n23  +  N°2 


N     =  Avogadro's  number 


=  0.6022 


cm^-  atoms 
b-mole 


M     =  Atomic  mass  of  indicated  nuclide,  [gm/mole] 

1  3 

p     =  Density  of  indicated  nuclide,  [gm/ciri  ] 

P.F.  =  Packing  factor  of  balls 
Volume  of  balls 


Volume  of  reactor 


0.61 


3 


V^  =   Volume  of  single  ball,  [cm  ] 

V'  =  Effective  volume  of  ball  (cell-averaged)  ,  [cm  "*  ] 

V  =  Volume  of  indicated  nuclide  in  ball,  [cm^] 

V  =  Effective  volume  of  nuclide  in  ball  (cell- 

averaged),  [cm^] 

2 

N  =  Atom  density  of  indicated  nuclide  in  ball, 

^      [b~'cm"M 

2 

N*  =  Effective  atom  density  of  indicated  nuclide 


X 


in  ball  [b  ^crn'M 


R,    =  Radius  of  ball,  [cm] 


Subscript  T,  C  or  U  indicates  Th^^^Oz,  graphite  or  U^^^Og, 
respectively. 

2 

Subscript  T,  C,  U  or  0  indicates  Th232,  graphite,  U233  or 

oxygen,  respectively. 
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D . 3 . 2   Material  Properties 

Calculation  of  initial  number  densities,  as  a  function  of 
U233  enrichment,  involved  considerations  of  the  volume  fractions 
of  the  particular  materials  within  a  fuel  cell.   Fuel  cells 
were  defined  as  the  volume  of  one  fuel  ball  divided  by  the 
I   packing  factor,  P.F.   Properties  of  the  major  substances  in 
the  fuel  balls  are  shown  in  Table  [D.2] 

The  makeup  of  the  fuel  balls  is  illustrated  in  Figure  [D.l], 
appropriately  dimensioned.   Packing  fractions  of  fueled  grains 
in  the  meat  are  listed  in  Table  [D.3] .   Number  densities  for 
all  materials  present  at  beginning-of-lif e  are  given,  and 
all  others  appearing  in  the  depletion  chains  are  initially 
zero. 
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TABLE  D . 2 
BALL  MATERIAL  PROPERTIES  [W3] 


Material 


Heavy  metal  oxides 
Graphite: 

grain  coating 

meat 

outer  shell 


X 


(gm/cm^ ) 


9.5 


M 


(gm/mole) 


264 


1.5 

12.0 

1.0 

12.0 

1.75 

12.0 
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3.00  cm 


BALL 

OUTER  SHELL 

MEAT 
GRAIN 

KERNEL 


FIGURE  D.l:   SECTIONAL  VIEW  OF  FUEL  BALL 
(not  to  scale) 
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TABLE  D . 3 
PACKING  FRACTIONS  OF  GRAINS  IN  MEAT  [W3] 


C/HM 

PACKING 

ATOM  RATIO 

FRACTION* 

110 

0.39073 

250 

0.168 

325 

0.127 

450 

0.091987 

.^   ,  .    ^    . .      (Number  of  grains) (Volume  per  grain) 
*Packing  fraction  =  -^^ —^ ^ 

Volume  of  Meat 
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D.3.3   Calculation  of  Number  Densities 
The  volume  of  a  fuel  element  is  given  by: 

Vg  -  |-7tR-^  =  113.1  cm^,  (D.l) 

so  the  effective  volume  (cell-averaged)  is: 

V^  =  pTjt;-  =  185.4  cm''.  (D.2) 

Number  densities  for  thorium,  oxygen,  uranium,  and  carbon  were 
determined  using  the  following  relationships: 

Phm^a 

PC^A 
N   = 
C     Mj^  .  (D.4 

%  =  -n^  (D.5) 

where  the  densities,  p,  are  "smeared"  throughout  the  ball. 
&  In  addition, 

1*       N        N 

C/™  =  N-f-N  =    N^  (°-^) 

T     U     HM 

N         N 
and      e  =     U     =   U   ,  (0.7)' 

N   +  N     N 
T     U     HM 

To  find  the  smeared  densities  for  heavy  metal  and  carbon: 

,  %M  U  ^Kernel  U^Grainlr  Meat  ,      ,^  on 


)2  /  \  Grain  j  \^    Meat  /\  Ball^ 
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Pc    ^Grainl-^    V 


Kernel  1    

Grain 


Grain    Meat 


V 


V 


"^  ^Meat  P    V. 


Grain 


.J 


Meat 


Meat 

Meat 
\^Ball| 


V 
"^  ^Outer  Shell  r  ~  V 


Meat 


Ball 


Ball 


(D.9) 


Since  N 


X 


N  (P.F.),  the  effective  number  densities  (cell- 


averaged)  can  be  summarized  as  follows 


N  ' 


N^' 


eo     I  HM 


M.,,,   \/V., A 


2  M 


HMO2/ 


Kernel 


V 


Ball 


N 


(P.F.) 


A 


M. 


HM 


M 


=  (1  -  e)P 


HM 


^^2  M^Mo 


N 


N, 


A 


N 


(P.F.) 


HM 


P^ 
Gram 


/V 


N. 


Kerne_l   (p.p.)  "A 


V 


Ball 


M 


HM 


V 


1  - 


Kernel 


V 


Grain 


/V 


Grain 
\  ^Ball 


(D.IO) 


(D.ll) 


■^  PMeat'l  - 


^Grain]  rMeat 
^Meatl  TBall 


I 


I   _   Meat 

^Outer  Shell      V,^  ,  t 

Ball/ 


(D.12) 


In  addition  to  the  above,  initial  number  densities  were 
used  as  input  for  silicon;   the  recvel  correction  (reciprocal 
velocity  [W3])  for  graphite  impurity  is  used  at  a  concentration 
of  approximately  2  ppb  carbon,  so  that  the  effective  number 
densities  (cell-averaged)  become: 


L 


N'      ,  =  (2.05E-9)N„' 
recvel  C 


N  '  ,  .    =8. 308E-5 
helium 


(D.13) 
(D.14) 
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^Meat 1 /^Grain 


N'.,.     =  (9.14E-3)L,  "-""  L,   """  (D.15) 

silicon  V--,  T  ,   V^,   , 

\  Ball / \  Meat  / 


Representative  values  for  N  '  and  N„,  '  are  given  in 

Table  [D.4],  while  values  for  N,'  ,  ,  N '  ,  N '  ,  N '     ., 

HM     O     C  '   recvel' 

N'.^.    ,  and  N'  ,.    are  listed  in  Table  [D.5]. 
silicon        helium 
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TABLE  D . 4 
INITIAL  CELL-AVERAGED  HEAVY  METAL  NUMBER  DENSITIES* 

e(%)    U233        Th232 

2.0   7.9408E-6  3.8910E-4 

4.5   1.7867E-5  3.7918E-4 

7.0   2.7793E-5  3.6925E-4 

2.0   3.4143E-6  1.6730E-4 

4.5   7.6821E-6  1.6303E-4 

7.0   1.1950E-5  1.5876E-4 

2.0   2.5810E-6  1.2647E-4 

4.5   5.8073E-6  1.2324E-4 

7.0   9.0336E-6  1.2002E-4 

2.0   1.8695E-6  9.1604E-5 

4.5   4.2063E-6  8.9267E-5 

7.0   6.5431E-6  8.6930E-5 


*  In  b   cm 


C/HM 

ATOM  RATIO 

110 

110 

110 

250 

250 

250 

325 

325 

325 

450 

450 

450 
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D. 4   Exposure  Times 

This  section  describes  the  calculation  of  exposure  times 
(hours)  used  in  PDQ-7  depletion  calculations  for  fuel  burnups 
up  to  120,000  MWD/MTHM.   Calculated  exposure  times  are  tabu- 
lated as  a  function  of  the  carbon-to-heavy  metal  atom  ratio 
in  Table  [D. 6]  . 

D.4.1   Nomenclature 

q' ' '  =  Average  core  power  density 

3 

=  5.0  watts/cm 

N'.,   =  Cell-averaged  heavy  metal  atom  density, 
HM  n        J 

[b   cm   ] 
t     =  Exposure  time  at  full  power,  [hours] 

Nj^    =  Avogadro's  number 

2 

=  0.6022  cm  -atoms/b-mole 

M,,,,   =  Atomic  mass  of  heavy  metal  nuclides 
HM 

-  2  32  gm/mole 

B     =  Fuel  burnup,  [MWD/MTHM] 

3. 
V^    =  Volume  of  relevant  reactor  portion,  [cm  J 
R 


-161- 


D . 4 . 2   Calculation  of  Exposure  Times 

Exposure  times,  t [hours],  were  calculated  to  correspond 

to  specific  values  of  fuel  burnup,  B  [MWD/MTHM] .   As  a  function 

of  the  cell-averaged  heavy  metal  number  density,  N'   [b   cm   ] , 

the  fuel  burnup  is  given  by: 

q'"N  t 
B  =   (2.4E-23)  ^^r-^ MWD/MTHM  ^    ^ 

HM  HM  Vi^.xu; 

so  that  the  exposure  time,  t  [hours] ,  can  be  expressed  as 

t  [hours]  =  (1. 8492E+3)B  [MWD/MTHM] N^^[b~"'"cm~-'-]   (D.17) 

where  N'   can  be  found  by  summing  Equations  (D.IO)  and  (D.ll). 
Representative  values  for  the  exposure  times  are  given  in 
Table  [D.6]  as  a  function  of  C/HM  atom  ratio. 
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TABLE  D . 6 
EXPOSURE  TIMES* 

C/HM  ATOM  RATIO 
B(MWD/MTHM)       110         250  325         450 

100      7.3421E+1   3.1568E+1  2.3864E+1  1.7285E+1 

1000      7.3421E+2   3.1568E+2  2.3864E+2  1.7285E+2 

2000      1.4684E+3   6.3136E+2  4.7728E+2  3.4570E+2 

3000      2.2026E+3   9.4704E+2  7.1592E+2  5.1855E+2 

6000      4.4053E+e   1.8941E+4  1.4318E+3  1.0371E+3 

10,000      7.3421E+3   3.1568E+3  2.3864E+3  1.7285E+3 

15,000      1.1013E+4   4.7352E+3  3.5796E+3  2.5927E+3 

20,000      1.4684E+4   6.3136E+3  4.7728E+3  3.4570E+3 

25,000      1.8355E+4   7.8920E+3  5.9660E+3  4.3212E+3 

30,000      2.2026E+4   9.4704E+3  7.1592E+3  5.1855E+3 

45,000      3.3040E+4   1.4206E+4  1.0739E+4  7.7783E+3 

60,000      4.4053E+4   1.8941E+4  1.4318E+4  1.0370E+4 

75,000      5.5966E+4   2.3676E+4  1.7898E+4  1.2964E+4 

90,000      6.6079E+4   2.8411E+4  2.1478E+4  1.5556E+4 

100,000      7.3421E+4   3.1568E+4  2.3864E+4  1.7285E+4 

120,000      8.8106E+4   3.7882E+4  2.8637E+4  2.0742E+4 

*Tiines  in  hours. 
For  comparison,  note  that  8766  hours  =  1  year. 


-163- 


D. 5   Depletion  Treatment 

D.5.1   Depletion  Chains 

Depletion  chains  for  the  46  nuclides  considered  are  shown 
in  Table [D. 7].   PDQ-7  solves  the  sequence  of  chains  in  order. 
Chain  1  through  Chain  19,  for  each  pass  through  the  depletion 
mesh.   In  Table  [D.7]  Chains  1  through  14  are  depleting. 
Underscoring  indicates  the  existence  of  a  fission  yield,  given 
in  Table  [D.9],  for  the  indicated  nuclide.   Reactions  are  in- 
dicated as  follows: 

(1)  :   (B") 

(2)  :   (n,Y) 

(3)  :   (n,Y)  plus  (6") 

D.5.2   Nuclide  Decay  Constants 

In  Table  [D.8]  are  tabulated  all  nuclide  decay  constants 
used  in  the  subject  PDQ-7  calculations. 

D . 5 . 3   Fission  Product  Yields 

Table  [D.9]  lists  the  yields  of  each  fission  product 
treated  for  six  fertile  and  fissile  nuclides. 

D . 5 . 4   Energy  Release  per  Fission 

The  energy  release  per  fission  for  all  heavy  metals  treated 
by  the  AMPX  code  are  given  in  Table  [D.IO]. 

D . 5 . 5  Chi  (Neutron  Yield)  Spectra 

As  a  function  of  the  carbon-to-heavy  metal  atom  ratio, 
the  groupwise  distributions  for  fission  neutrons,  X/  are  given 
in  Table  [D. 11] . 
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TABLE  D.7  (page  1  of  4 ) 
DEPLETION  CHAINS 


CHAIN  1:   U233  Chain 


Pa233 

(2)       (2)       (2) 
U2  33  >-U234 >-U23  5 ^U236 


CHAIN  2:   Xel35  Chain 

1135 

(1)    \ 

Xel35 

(1)    \ 

Csl35 

CHAIN  3:   Sml49  Chain  1 


Ndl47 


(1)    \  (2)         (2) 

Pml47 >  Pml48m — ^  Pml49 


(3) 


(2) 


Sml52 
(3' 


(2) 


(2) 


Sml49 >  SmlSO >-Sml51 


(2)                    (2) 
Eul53 ^Eul54 >■  Eul55 
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TABLE  D.7  (page  2  of  4 ) 


CHAIN  4:   Sml49  Chain  2 


Ndl47 


(1) 


Pm 


(2)         (2) 

4  7 >  Pml4  8  >  Pml49 

(3)^ 


(2)         (2) 
Sml49 >  SmlSO  >  SmlSl 


(2) 
— ^Sml52 


(3) 


(2)  (2) 

-Eul53 >Eul54 >Eul55 


CHAIN    5:       Kr83    Chain 

(2) 

Kr82 >Kr83 


CHAIN  6:   Rhl03  Chain 


Rhl0  3 


(1) 


CHAIN  7:   Rhl05  Chain 


Rhl05 


(1) 


CHAIN  8:   Agl09  Chain 


Agio  9 


(1) 
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TABLE  D.7  (page  3  of  4  ) 


CHAIN  9:   Xel31  Chain 


XelSl 


(1) 


CHAIN  10:   Fission  Product  Chain  (lumped) 

(2) 

FPl >■  FP2 


CHAIN  11:   Ndl43  Chain 


Ndl43 


(1) 


CHAIN  12:   Xel33  Chain 


Xel33 


(1) 


(2) 
Csl33  3-CS134 


CHAIN  13:   Ndl45  Chain 


Ndl45 


(1) 


CHAIN  14:   Pu239  Chain 


U238 
(3) 


^Np239 

^■^^  I  (2)         (2)         (2) 

Pu239 ^  Pu24  0 ^  Pu241 :^Pu24  2 

(3r 


'Am243 
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TABLE  D.7  (page  4  of  4) 
CHAIN  15:   Helium  Chain  (non-depleting) 

Helium 
CHAIN  16:   Carbon  Chain  (non-depleting) 

Carbon 
CHAIN  17:   Oxygen  Chain  (non-depleting) 

Oxygen 

V   CHAIN  18:   Recvel  Chain  (non-depleting) 

I 

! 

(  Recvel 

.    CHAIN  19:   Silicon  Chain  (non-depleting) 
Silicon 
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TABLE  D.8  (page  1  of  2) 
NUCLIDE  DECAY  CONSTANTS 


CHAIN 

I# 

REFERENCE 

NUCLIDE 

*X, [sec   J 

# 

NUMBER 

1 

23 

23290 

Th232 

1.58E-18 

1 

42 

23391 

Pa233 

2.93E-7 

1 

24 

23392 

U233 

1.358E-13 

1  . 

25 

23492 

U234 

8.856E-14 

1 

25 

23592 

U235 

3.08E-17 

1 

18 

23692 

U236 

9.189E-16 

2 

4 

13553 

1135 

2.87E-5 

2 

8 

13554 

Xel35 

2.09E-5 

2 

11 

13555 

Csl35 

1.098E-14 

3,4 

3 

14760 

Ndl47 

7.32E-7 

3,4 

12 

14761 

Pml47 

8.29E-9 

3 

13 

14800 

Pml48m 

1.98E-7 

4 

14 

14861 

Pml4  8 

1.49E-6 

3,4 

28 

14961 

Pml4  9 

3,63E-6 

3,4 

15 

14962 

Sml49 

5.491E-23 

3,4 

38 

15062 

Sml50 

Stable 

3,4 

16 

15162 

Sml51 

2.75E-10 

3,4 

39 

15262 

Sml52 

Stable 

3,4 

17 

15363 

Eul53 

Stable 

3,4 

40 

15463 

Eul54 

1.37E-9 

3,4 

41 

155G3 

Eul55 

1.16E- 8 

5 

29 

8236 

Kr8  2 

Stable 

5 

30 

8336 

Kr83 

Stable 

6 

31 

10345 

Rhl03 

Stable 

7 

32 

10545 

Rhl05 

5.36E-5 

8  " 

33 

10947 

Agl09 

Stable 
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TABLE  D.8  (page  2  of  2) 


9 

34 

13154 

Xel31 

Stable 

10 

22 

10001 

FPl 

Stable 

10 

21 

10002 

FP2 

Stable 

11 

36 

14360 

Ndl43 

Stable 

12 

35 

13354 

Xel33 

1.52E-6 

12 

9 

13355 

Csl33 

Stable 

12 

10 

13455 

Csl34 

9.55E-9 

13 

37 

14560 

Ndl45 

Stable 

14 

43 

23892 

U238 

4.87E-18 

14 

27 

23993 

Np239 

3.41E-6 

14 

44 

23994 

Pu249 

9.016E-13 

14 

45 

24094 

Pu240 

3.249E-12 

14 

46 

24194 

Pu241 

1.68E-9 

14 

19 

24294 

Pu242 

5.795E-14 

14 

20 

24395 

Ain243 

2.871E-12 

15 

1 

402 

Helium 

Stable 

16 

2 

1206 

Carbon 

Stable 

17 

5 

1608 

Oxygen 

Stable 

18 

6 

1200 

Recvel 

Stable 

19 

7 

1428 

Silicon 

Stable 

*  For  stable  nuclides,  X    =    O.lE-20  sec    was  used 
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TABLE  D.9 
FISSION  PRODUCT  YIELDS 

FISSION  I#  FISSION  SOURCE 

PRODUCT       Th232     U233      U235      U238      Pu239     Pu241 


T   1135 

4 

0.0539 

0.0562 

0.0617 

0.0578 

0.0693 

0.0626 

Xel35 

8 

0.0021 

0.0022 

0.0024 

0.0022 

0.0027 

0.0024 

Pml4  9 

28 

0.00945 

0.0077 

0.0113 

0.021 

0.013 

0.012 

Ndl43 

36 

0.08 

0.0591 

0.0603 

0.48 

0.046 

0.056 

FPl 

22 

0.301 

0.309 

0.298 

0.273 

0.324 

0.324 

FP2 

21 

1.137 

1.119 

1.26 

1.426 

1.456 

1.456 

Kr82 

29 

0.01 

0.0069 

0.0028 

0.0005 

0.0017 

0.0007 

Kr83 

30 

0.0199 

0.0114 

0.00544 

0.0009 

0.0029 

0.001 

Rhl03 

31 

0.0016 

0.016 

0.03 

0.062 

0.056 

0.062 

RhlOS 

32 

0.0005 

0.005 

0.009 

0.037 

0.055 

0.058 

Agl09 

33 

0.00055 

0.0004 

0.0003 

0.003 

0.0155 

0.031 

Xel31 

34 

0.017 

0.0341 

0.0293 

0.023 

0.038 

0.034 

Xel33 

35 

0.037 

0.0588 

0.659 

0.052 

0.069 

0.06 

Ndl45 

37 

0.066 

0.0338 

0.0398 

0.039 

0.031 

0.036 

Eul53 

17 

0.00038 

0.0015 

0.00169 

0.0038 

0.0034 

0.002 

Eul55 

41 

0.00009 

0.0003 

0.00033 

0.0013 

0.0016 

0.001 

Ndl47 

3 

0.038 

0.0193 

0.0236 

0.028 

0.0205 

0.022 

SmlSl 

16 

0.002 

0.0035 

0.0044 

0.009 

0.008 

0.005 

Sml52 

39 

0.00085 

0.0022 

0.00281 

0.006 

0.006 

0.003 
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TABLE  D.IO 
ENERGY  RELEASE  PER  FISSION 


FISSIONABLE 

< 

NUCLIDE 

[watt-sec/fission 

Th232 

3.18 

X 

10"^  ^ 

Pa233 

3.19 

X 

10"^  ^ 

U233 

3.19 

X 

10"^  ^ 

U234 

3.19 

X 

10"^  - 

U235 

3.21 

X 

10"^^ 

U236 

3.26 

X 

10"^  ^ 

U238 

3.31 

X 

10"^  ^ 

Np239 

3.27 

X 

io~-  - 

Pu239 

3.34 

X 

10~^  ^ 

Pu240 

3.36 

X 

10~^  ^ 

Pu241 

3.37 

X 

10"'  ^ 

Pu242 

3.37 

X 

io~-  - 

Ain243 

3.39 

X 

10"^  ' 
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TABLE  D.ll 

Cni  (NEUTRON  YIELD)  SPECTRA 

C/HM  ATOM  RATIO 
110           250  325  450 

Xj^     9.66727E-1    9.66727E-1  9.66727E-1    9.66727E-1 

X2     3.32743E-2    3.32745E-2  3.32746E-2    3.32743E-2 

X3     4.90742E-9    5.67340E-9  5.91979E-9    5.17704E-9 

X4     4.98079E-14   5.75821E-14  6.00830E-14   5.25449E-14 
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D . 5 . 6   Neutron  Flux 

Table  [D.12]  tabulates  the  groupwise  relative  neutron 
flux  generated  for  each  carbon-to-heavy  metal  atom  ratio  by 
the  AMPX  code.  Group  parameters  are  given  in  Table  [D.13]. 
These  flux  values  were  used  to  collapse  the  four-group  con- 
stants given  in  Appendix  C  to  one-group,  for  use  in  the 
analysis  of  Chapter  2. 
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TABLE  D.12 
GROUPWISE  RELATIVE  NEUTRON  FLUX* 

C/HM  ATOM  RATIO 
FLUX    110     250     325     450 

(f)-j^  1.000  1.000  1.000  1.000 

(|)2  1.534  1.533  1.533  1.534 

(J)3  1.025  1.194  1.236  1.282 

(})^  0.920  2.305  2.939  3.874 

*  Normalized  by  using  ^    /cj)-,  for  each  C/HM 
atom  ratio. 
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TABLE  D.13 
ENERGY  GROUPS 


Group  1  Group  2  Group  3  Group  4 

Upper  Energy  (eV)    1.4918E+7  1.8316E+5  5.8295E+2  1.8554E+0 

Lower  Energy  (eV)    1.8316E+5  5.8295E+2  1.8554E+0  5.00E-3 

Lethargy  width,  Au   4.400  5.750  5.750  5.916 
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